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Abstract 



The production of protons, anti-protons, neutrons, deuterons and tritons in minimum 
bias p+C interactions is studied using a sample of 385 734 inelastic events obtained 
with the NA49 detector at the CERN SPS at 158 GeV/c beam momentum. The data 
cover a phase space area ranging from to 1.9 GeV/c in transverse momentum and in 
Feynman x from -0.80 to 0.95 for protons, from -0.2 to 0.4 for anti-protons and from 
0.2 to 0.95 for neutrons. Existing data in the far backward hemisphere are used to 
extend the coverage for protons and light nuclear fragments into the region of 
intranuclear cascading. The use of corresponding data sets obtained in hadron-proton 
collisions with the same detector allows for the detailed analysis and 
model-independent separation of the three principle components of hadronization in 
p+C interactions, namely projectile fragmentation, target fragmentation of participant 
nucleons and intranuclear cascading. 



1 Introduction 



Baryon and light ion production in proton-nucleus collisions has in the past drawn con- 
siderable interest, resulting in an impressive amount of data from a variety of experiments. This 
interest concentrated in forward direction on the evident transfer of baryon number towards the 
central region, known under the misleading label of "stopping", and in the far backward region 
on the fact that the laboratory momentum distributions of baryons and light fragments reach 
far beyond the limits expected from the nuclear binding energy alone. A general experimental 
study covering the complete phase space from the limit of projectile diffraction to the detailed 
study of nuclear effects in the target frame is, however, still missing. More recently, renewed 
interest has been created by the necessity of providing precision reference data for the control 
of systematic effects in neutrino physics. 

In addition to and beyond the motivations mentioned above, the present study is part of 
a very general survey of elementary and nuclear interactions at the CERN SPS using the NA49 
detector, aiming at a straight-forward connection between the different reactions in a purely 
experiment-based way. After a detailed inspection of pion [1], kaon [|2l and baryon [3 J produc- 
tion in p-i-p interactions, a similar in-depth approach is being carried out for p-i-C collisions. This 
has led to the recent publication of two papers concerning pion production flUO and this aim is 
here being extended to baryons and light ions. 

The use of the light, iso-scalar Carbon nucleus is to be regarded as a first step towards 
the study of proton collisions with heavy nuclei using data with controlled centrality available 
from NA49. It allows the control of the transition from elementary to nuclear interactions for a 
small number of intra-nuclear collisions, thus providing an important link between elementary 
and multiple hadronic reactions. It also allows for the clean-cut separation of the three basic 
components of hadronization in p-i-A collisions, namely projectile fragmentation, fragmenta- 
tion of the target nucleons hit by the projectile, and intra-nuclear cascading. The detailed study 
of the superposition of these components in a model-independent way is the main aim of this 
paper. For this end the possibility of defining net proton densities by measuring anti-protons 
and thereby getting access to the yield of pair produced baryons, will be essential. As the ac- 
ceptance of the NA49 detector does not cover the far backward region, the combination of the 
NA49 results with measurements from other experiments dedicated to this phase space area 
is mandatory. A survey of the s-dependence of backward hadron production in p-i-C collisions 
has therefore been carried out and is published in an accompanying paper This allows the 
extension of the NA49 data set to full phase space. 

As the extraction of hadronic cross sections has been described in detail in the preced- 
ing publications IIH3, the present paper will concentrate on those aspects which are specific 
to baryon and light ion production, particularly in the exploitation of the NA49 acceptance 
into the backward hemisphere. After a short comment on existing double differential data in the 
SPS energy range in Sect. [21 a few experimental details will be given in Sect.[3]together with the 
binning scheme adopted for protons, anti-protons and neutrons. Section|4]will present a compre- 
hensive description of particle identification in the backward hemisphere which is an important 
new ingredient of the optimized use of the NA49 detector in particular for the asymmetric p-i-A 
collisions. Section [5] deals with the extraction of the inclusive cross sections and with the ap- 
plied corrections. Section [6] contains the data tables and plots of the invariant cross sections as 
well as some particle ratios and a comparison to the few available double differential yields at 
SPS energy for comparison. Section|7]describes the use of the extensive complementary data set 
from the Fermilab experiment [7J for the data extension into the far backward direction together 
with an interpolation scheme allowing for the first time the complete inspection of the produc- 
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tion phase space for protons in the range -2 < xp < +0.95. This combined study is extended 
to deuterons and tritons in Sect. [H Baryon ratios are shown in Sect. |9l px integrated quantities 
are given in Sect. [10] both for minimum bias trigger conditions and for the dependence on the 
number of measured grey protons. In addition, the measured integrated neutron yields are 
presented in Sect. [10] together with a comparison to other integrated data in the SPS energy 
range. Section [TT] contains a detailed discussion of the two-component mechanism of baryon 
and baryon pair production, thus covering the first two components of the hadronization process 
defined above. This Section contains experimental results from both p+p and p+C interactions 
including a comment on resonance decay and a comparison to a recent microscopic simulation 
code. Section [12] gives a detailed discussion of anti-proton production including the application 
of the two-component mechanism introduced in Sect.[lT]and a study of the px dependence. The 
discussion of pt integrated proton and net proton yields is presented in Sect. [T3]followed by the 
exploitation of double differential proton and net proton cross sections in Sect. [14] The paper is 
closed by a summary of conclusions in Sect. [TS] 

2 The Experimental Situation 

As already pointed out for pions in [4J there are only two sets of double differential 
inclusive data, for identified baryons and light fragments in p+C collisions in the SPS energy 
range. The differential inclusive cross sections are presented in this paper as: 



dx pdp^ ' 

with xp = 2pl/^ defined in the nucleon-nucleon cms. A first data set [T'.'Sl covers the far 
backward direction for protons and light ions at five fixed laboratory angles between 70 and 
160 degrees for total lab momenta between 0.4 and 1.4 GeV/c at a projectile momentum of 
400 GeV/c. A second set [9] has been obtained in forward direction for 0.3 < xp < 0.88 and 
0.15 < Pt < 0.5 with 100 GeV/c beam momentum. The respective phase space coverage in xp 
and Pt is shown in Fig.[T^,b for protons and anti-protons, respectively, with a superposition of 
the NA49 coverage for protons. This coverage is presented in more detail in Fig. [it and for the 
anti-protons in Fig.[T]i. 

With the NA49 data covering lab angles of up to 40 degrees the combination with l[71[8l 
into a consistent data set becomes possible. This allows for the first time the complete scrutiny of 
the proton phase space in the range -2 < xp < +0.95, with only minor inter- and extrapolation. 

3 Experimental information and binning scheme 

As a detailed description of the NA49 detector and the extraction of inclusive cross 
sections has been given in [[Tl-[4[[T0l. only some basic informations are repeated here for conve- 
nience. 

3.1 Target, grey proton detection, trigger cross section and event sample 

The NA49 experiment is using a secondary proton beam of 158 GeV/c momentum at the 
CERN SPS. A graphite target of 1.5% interaction length is placed inside a grey proton detector 
which covers a range from 45 to 315 degrees in polar angle with a granularity of 256 readout 
pads placed on the inner surface of a cylindrical proportional counter. An interaction trigger is 
defined by a small scintillator 380 cm downstream of the target in anti-coincidence with the 
beam. This yields a trigger cross section of 210.1 mb corresponding to 91% of the measured 
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Figure 1 : Phase space coverage of existing data: a) p data from f7] (full lines) and BUl. Here with 
the shaded area is shown the NA49 acceptance range; b) p data from [9l; c) p data from NA49; 
and d) p data from NA49 



inelastic cross section of 226.3 mb. This is in good agreement with the average of 225.8 mb 
obtained from a number of previous measurements [4J. A total sample of 385.7k events has 
been obtained after fiducial cuts on the beam emittance and on the longitudinal vertex position. 

3.2 Acceptance coverage, binning and statistical errors 

The NA49 detector [[TOll covers a range of polar laboratory angles between ±45 degrees 
with a set of four Time Projection Chambers combining tracking and particle identification, two 
of the TPC's being placed inside superconducting magnets. While for anti-protons the accessi- 
ble range in xp and px is essentially defined by the limited event statistics, it has been possible 
to completely exploit the available range of polar angle for protons. The corresponding binning 
schemes are shown in Fig.|2]in the cms variables xp and pr- 

A rough indication of the effective statistical errors is given by the shading of the bins. 
Neutrons have been detected in a forward hadronic calorimeter [3] in combination with propor- 
tional chambers vetoing charged hadrons. Due to the limited resolution in transverse momentum 
only pt integrated information in 8 bins in xp (Fig. |2]:) could be obtained, after unfolding of 
the energy resolution. This coverage is identical to the one in p-i-p interactions ^ and allows 
for direct yield comparison. 
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Figure 2: Binning scheme inxp and px together with information on the statistical errors for a) 
protons, b) anti-protons and c) neutrons 



4 Particle Identification 

Due to the forward-backward asymmetry of p-i-A interactions, the study of the backward 
hemisphere is of major interest for the understanding of target fragmentation and intranuclear 
cascading. Particle identification at negative Xi;' is therefore mandatory; it has to rely for the 
NA49 detector on the measurement of ionization energy loss in the TPC system. This method 
has been developed and described in detail for mesons and baryons in p-i-p collisions in [[l]-[3l 
for xp > 0. A substantial effort has been invested for the present study in its extension to 
the far backward direction down to the acceptance limit in Xi;' imposed by the NA49 detector 
configuration. With decreasing x p the baryonic lab momentum decreases below the region of 
minimum ionization where the ionization energy loss increases like and thereby succes- 
sively crosses the deposits from kaons, pions and electrons. This is shown in Fig.[3]for the the 
momentum dependence of the parametrization of the mean truncated energy loss used in this 
analysis. 




Figure 3: Parametrization of the mean truncated energy loss as a function of total lab momen- 
tum piab for electrons, pions, kaons and protons. The situation for deuterons and tritons is also 
indicated 
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In terms of xp and pr, this cross-over pattern reflects into lines of equal energy loss in 
the Xp-PT plane as shown in Fig. HI 




Figure 4: Lines of equal energy loss for protons and kaons (p-K), protons and pions (p-vr) and 
protons and electrons (p-e) as functions of xp and px, together with the acceptance limit of the 
NA49 detector 

The region above the line p-K allows for the standard multi-parameter fits of the trun- 
cated energy loss distributions as described in the preceding publications [[I]-31|. The approxi- 
mately triangular region below the line p-e permits the direct extraction of baryon yields par- 
tially even without fitting. This is exemplified in Fig. [5] for two bins at xp = -0.5 and -0.6 and 
small pt. 
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Figure 5: Truncated energy loss distributions for a) positives and b) negatives at xp = -0.5, 
= 0.1 GeV/c; and for c) positives and d) negatives at xp = -0.6, pr = 0.15 GeV/c 



5 



It is interesting to note that also the light ions deuteron and triton are here well separable 
practically without background. For anti-protons, a direct measurement of the p/p ratio becomes 
feasible down to values below 10^'^ in this region. 

In order to extract proton yields from the energy loss distributions in the intermediate 
region between the lines p-e and p-K, Fig HI the particle ratios in each studied xp/pr bin are 
of prime importance. If the ratios p/vr, k/vr and e/vr are known, proton yields may be obtained 
from the total number of tracks even in those bins where the proton energy loss equals the one 
from electrons, pions or kaons. A two-dimensional interpolation of the measured particle ratios 
over the full accessible phase space has therefore been established. These ratios are obtained 
without problem in the regions below the line p-e and above the line p-K (Fig. H]) as well as in 
most intermediate bins where a sufficient separation in dE/dx of the different particle species 
is present. Near the cross-over bins the measured ratios show a sharp increase of the effective 
statistical fluctuations, an increase which has been described in the discussion of the error matrix 
involved with the multi-dimensional fitting procedure in [2]. As this effect is of statistical and 
not of systematic origin, an interpolation through the critical regions in x p and pt is applicable. 

It should be stressed here that the obtained particle ratios are non-physical in the sense 
that they use different phase space regions for each particle mass. For each xplpr bin the nec- 
essary transformation to total lab momentum is performed using the proton mass for each track. 
This means that electrons, pions and kaons from different effective xp values enter into the 
proton bin, with an asymmetry that increases with decreasing xp and pt- This is quantified in 
Fig. |6] where the effective mean xp for electrons, pions and kaons is shown as a function of 
proton x^p for two values of pt- 
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Figure 6: xp for electrons, pions and kaons as a function of proton x^p for pp = Q.l GeV/c (upper 
lines) and Pt= 1-3 GeV/c (lower lines) 



The lighter particles at small pr are thus effectively collected from the neighbourhood 
of Xp = with decreasing baryonic xp. For anti -protons this purely kinematic effect is un- 
favourable for fitting as the effective p/7r^ ratios quickly decrease below the percent level at low 
Pt, whereas the p/7r+ ratios stay always above about 10%, increasing rapidly with px due to the 
rather flat number distribution dn/dxp. Proton and anti-proton extraction are therefore regarded 
separately in the following Sects. I4n and l4.2[ respectively. 
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4.1 



Proton extraction 



4.1.1 e'^/n^ ratio 

The crossing of the proton dE/dx through the practically constant electron energy loss 
at piab ~ 1 GeV/c is the least critical effect as the momentum dependence of the proton dE/dx 
is a steep function of lab momentum in this piab range and as the e/vr ratio quickly decreases with 
increasing pt, reaching the 1% level already at pr > 0.3 GeV/c. The e+/7r+ ratio is shown in 
Fig- Has a function of xp for four values of px together with the two-dimensional interpolation 
used. 
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Figure 7: e+/7r+ ratio as a function of xp for a) pt = 0.1 GeV/c, b) px = 0.2 GeV/c, c) 
Pt = 0.3 GeV/c and d) px = 0.4 GeV/c. The cross-over with the proton energy loss is indi- 
cated by vertical dotted lines. The full lines give the two-dimensional interpolation established 
from the data 



Vertical lines indicate the position of the dE/dx cross-over for each pt value and evi- 
dently the ratios may be well interpolated through the small affected xp regions. 

4.1.2 p/n^ ratio 

Fitted p/7r+ ratios are presented in Fig. [8] for four x p values together with their two- 
dimensional interpolation as a function of pp. If the fit results yield stable pp dependences 
within their statistical errors in the uncritical regions at Xi? = and -0.6, the intermediate xp 
values at -0.2 and -0.4 show some additional fluctuation in the cross-over regions indicated by 
the hatched areas which combine the p-vr and p-K ambiguities. 

Evidently the data interpolation describes the ratio properly through the ambiguous pp 

areas. 

A complete picture over the full available backward phase space is given in Fig. |9] where 
the fitted p/7r+ ratios are shown as functions of xp for different pp values together with the 
interpolations (full lines). The ratios at successive pp values are shifted by a factor of 2 for 
clarity of presentation. 

The complete situation for the data interpolation is finally presented in Fig. [TO] with fixed 
vertical scale as a function ofxp at different pp values. Here the acceptance limit of the NA49 
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Figure 8: p/Tr"*" ratios as a function of pt for a) xp = 0, b) xp = -0.2, c) xp = -0.4 and d) = - 
0.6. The full lines present the two-dimensional data interpolation, the hatched areas between the 
vertical lines the regions affected by the p-n and p-K ambiguities 
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Figure 9: p/n^ ratios as functions of a;^ for fixed values of px [GeV/c]. Full lines: data interpo- 
lation. The ratios at successive px values are shifted by a factor of 2 for clarity of presentation 

detector is given as the broken line together with the region of p-vr and p-K ambiguity as hatched 
area. 

This plot again clarifies the way in which the critical cross-over areas may be bridged 
by two-dimensional interpolation. 



4.1.3 ratio 

A situation quite similar to the p/n^ ratio exists for the K+/7r+ ratio. Again, there are 
regions of ambiguity against protons and pions, but the influence of eventual systematic devi- 
ations on the extraction of protons is small as the K+/7r+ ratios are smaller than the p/Tr"*" ratio 
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Figure 10: Interpolated p/7r+ ratios as functions of xp for fixed values of pr [GeV/c]. Broken 
line: NA49 acceptance limit. Hatched area: region of p-n and p-K ambiguity 



by factors between 3 and 10. Fig. [TT] shows K+/7r+ ratios as functions of for four xj? values, 
where the lowest and highest xj? at -0.6 and allow for unambiguous fits over the full xp range, 
whereas the xp values at -0.2 and -0.4 suffer from p-K and K-tt ambiguities in the hatched 
areas of px, with resulting increased statistical fluctuations. The two-dimensional interpolation 
is superimposed as full lines. 
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Figure 11: K+/7r+ ratios as a function of px for four values of xp, a) xp = 0, b) xp = -0.2, c) 
xp = -0.4 and d) xp = -0.6. The regions of p-K and K-vr ambiguities are indicated as hatched 
areas in panels b) and c). The full lines represent the two-dimensional interpolation 

All fitted values of K+/7r+ are plotted in Fig. [121 as a function of xp for fixed px- As 
in Fig. |9] the ratios at successive pp values are shifted by 3 in order to sufficiently separate the 
measurements. 
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Figure 12: K+/7r+ ratios as functions of for fixed values of pt [GeV/c]. Full lines: data inter- 
polation. The ratios at successive pt values are shifted by a factor of 3 for clarity of presentation 

Fig.[T3]presents the overview of the interpolated K+Ztt"*" ratios at fixed vertical scale as a 
function oixp for fixed values of px [GeV/c]. The broken line represents the acceptance limits 
and the hatched area the region of p-K and K-tt ambiguity. 
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Figure 13: Interpolated K+/7r+ ratios as functions of xp for fixed values of px- Broken line: 
NA49 acceptance limit. Hatched area: region of K-vr and p-K ambiguity 



4.1.4 Proton extraction in the far forward region 

Due to the gap between the TPC detectors imposed by the operation with heavy ion 
beams ifTOl . charged particles progressively leave the TPC acceptance region at low pt for 
xp > 0.55. Here, tracking is achieved by the combination of a small "gap" TPC (GTPC) in 
conjunction with two forward proportional chambers (VPC). The performance of this detector 
combination is described in detail in [3]. In the absence of particle identification in this area one 
has to rely on external information concerning the combined fraction of K+ and 7r+ in the total 
charged particle yield. Several considerations help to establish reference values for the (K+ -i- 
7r+)/p ratios: 
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- The (K+ + 7r+)/p ratios decrease very rapidly with increasing xp at all pr, from about 
10% at xp = 0.6 to less than 1% at = 0.9. Possible deviations from the used external 
reference data therefore introduce only small systematic effects in the extracted proton 
yield. 

- Existing data may be used to come to a consistent estimation of the particle ratio. Direct 
measurements from Barton et al. flU in p+C interactions cover the region from xp = 0.2 
to 0.8 for = 0.3 and 0.5 GeV/c. Although the published invariant cross sections show 
sizeable deviations from the NA49 results, see Sect. 16.51 the particle ratios of the two 
experiments compare well. 

- The ratios also comply with measurements in p+p collisions, both from NA49 |[l]-[3l and 
from Brenner et al. [fTTTl at 100 and 175 GeV/c beam momentum. 

An overview of the experimental situation is given in Fig. [14] which shows the available 
measurements of the (K+ + 7r+)/p ratio as a function of xp for eight values of px between 0.1 
and 1 .3 GeV/c. Here the full lines and full circles represent the NA49 measurements in p+p and 
p+C collisions, respectively, and the open circles and squares show the corresponding ratios 
from [H and [llj. 
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Figure 14: Measured (K"*" + 7r"'")/p ratios as a function of xp for eight values of px between 
0.1 and 1.3 GeV/c. The full lines give the results from NA49 in p+p collisions, the full and 
open circles the ones from NA49 and [9], respectively. Open squares: results from ifTTTl in p+p 
interactions 



An impressive consistency between the different experimental results within their sta- 
tistical uncertainties is apparent. This allows for the safe extrapolation of the ratios into the 
non-measured region in p+C interactions. The slight deviations towards xp = 0.3 are due to the 
fact that the reference ratios are obtained from invariant cross sections, whereas the values of in- 
terest for proton extraction are rather density ratios using the proton mass for the transformation 
to lab momentum, see Fig. [6] above. At a;i7' > 0.3 the two definitions converge rapidly. Also at 
Pt > 1.1 GeV/c systematic differences on the level of one standard deviation might be present. 
In this region however, the full TPC information is available at all xp values, see Fig. [15] which 
shows the interpolated (K+ + 7r+)/p ratios as a function of px for different values of xp. 

The general data consistency on a ~10% level apparent in Fig. [14] implies that possible 
systematic effects on the extracted proton cross sections in the region between the dotted and 



11 




[GeV/c] 

Figure 15: Interpolated (K"*" + 7r+)/p ratios as a function of px for different values of xp be- 
tween 0.3 and 0.9, full lines. Broken line: border between available TPC information and the 
GTPC/VPC combination (Tracking only). Dotted line: acceptance limit of the NA49 detector 



broken lines in Fig. [15] are on the percent level or below. 

4.2 Anti-proton extraction 

As stated above the extension of the determination of anti-proton yields into the back- 
ward hemisphere suffers from the fact that the p/vr and p/K^ ratios come down rapidly with 
decreasing xp- This is shown by the energy loss distributions of two typical bins in xp and pt 
in Fig.[T6l 
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Figure 16: dE/dx distributions for negative particles a) x^? = -0.1, = 0-1 and b) xp = -0.15, 
Pt = 0.3 

This effect is largely due to the asymmetry between the effective xp for light particles 
and anti-protons due to the transformation to the lab momentum using proton mass, see Fig.[6j 
Thus at x^p = -0.2 pions are sampled close to maximum yield whereas the anti-proton cross 
section is steeply decreasing. 

If the extraction of pion yields therefore presents no problem in this phase space region, 
the fits of kaon and anti-proton densities become strongly correlated with sizeable uncertainties 
in their relative position on the energy loss scale. The combined (K^ -i- p) yields however stay 
well defined with respect to the pions. This is shown by the fitted (K~ -i- ratios of Fig.fTTl 

In order to resolve the K -p ambiguity, the high statistics data on p production in p-i-p 
interactions [O have been invoked. In this symmetric configuration, the measured p cross sec- 
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Figure 17: (K^ + p)/7r^ ratios as a function of x^? for fixed values of pt- Full lines: two- 
dimensional interpolation of the fitted ratios 

tions may be reflected into the backward hemisphere and thereby the correlation between the 
relative shifts S-p and 5k- the dE/dx variable may be obtained imposing for each combination 
the given reflected cross section. As shown in Fig.[T8]for xp = -0.1 this correlation is precisely 
defined experimentally over the full px range, also yielding stable values of the p/ir^ ratio. 

^ 0.03 
0.02 
0.01 
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^ 

Figure 18: Correlation between the relative displacements 6p and 6^- in p-i-p collisions at x^? = - 
0.1 for fixed values of Pt, imposing the forward-backward symmetry of cross section in this 
interaction. The lines are given to guide the eye 

Using the same correlation for the p-i-C data, effective p/vr" ratios are obtained. The 
observed stability of these ratios over the full range of the correlations is a strong test of the 
validity of the method. 

Fig. (represents the obtained p/vr^ ratios as a function of pt for xp = -0.05, -0.1 and 
-0.15, together with the directly fitted values at xj? = 0. Full lines: two-dimensional interpolation 
of the ratios. 

A complete picture of the p/vr" ratios used in this analysis is given in Fig. [20] as a 
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Figure 19: p/vr ratios as a function of pr for a) xp = 0, b) xp = -0.05, c) xp = -0.1 and d) 
= -0.2 



function of which combines the directly fitted ratios in the forward hemisphere with the 
ones obtained using the reflection method described above in the backward hemisphere. 
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Figure 20: p/vr" ratios as a function of s^;^ for fixed values of -pp. The values at Xi? = -0.2 are 
extrapolations using the broken lines 



Note that the values "dX. xp = -0.2 are obtained by extrapolation following the broken 
lines. Note also that the applied method allows the extraction of the ratios in the percent and 
sub-percent region. 

5 Evaluation of invariant cross sections and corrections 

The invariant cross section, 

f{xF,PT) = E{xp,PT)--^{xp,PT) (2) 
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is experimentally determined by the measured quantity ^ 



fmcUxF,PT,Ap ) = E{XF,PT,AP )■— ^ , (3) 



o'trig An{xF,PT , Ap^ 

where Ap^ is the finite phase space element defined by the bin width. 

As described in [IJ several steps of normalization and correction are necessary in order 
to make fmeas{xF,PT, Ap^) approach f{xF,PT)- The determination of the trigger cross section 
and its deviation from the total inelastic p+C cross section have been discussed in [|4||. The 
following corrections for baryons have been applied and will be discussed below: 

- treatment of the empty target contribution 

- effect of the interaction trigger 

- feed-down from weak decays of strange particles 

- re-interaction in the target volume 

- absorption in the detector material 

- effects of final bin width 



5.1 Empty target contribution 

This correction has been determined experimentally using the available empty target 
data sample, as described in [4]. The resulting correction is essentially determined by the dif- 
ferent amounts of empty events in full and empty target condition. It is within errors pt inde- 
pendent and equal for protons and anti-protons. It increases from about 2% in in the far forward 
direction to about 7% in the most backward region as shown in Fig. [2lJ 




Figure 21: Empty target correction as a function of xf- The full line shows the chosen interpo- 
lation 



5.2 Effect of the interaction trigger 

Due to the high trigger efficiency of 93% this correction is small compared to p-i-p 
interactions [IJ. It has been determined experimentally by increasing the diameter of the trigger 
counter using the accumulated data. Within its statistical uncertainty it is independent of pt and 
similar for protons and anti-protons. Its xf dependence as shown in Fig. |22]is following the 
expected trend ^ where the fast decrease in forward direction as compared to p-i-p collisions is 
due to the lower particle yields at high xf- 
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Figure 22: Trigger bias correction as a function of xp for a) protons and b) anti-protons. The 
chosen interpolation is given by the full lines 

5.3 Feed-down correction 

The hyperon cross sections relative to p+p collisions established in llH for > have 
been used. These cross section ratios approach atxp < the expected factor of 1.6 correspond- 
ing to the number of intranuclear projectile collisions, in account of the fact that for A and A 
there is no isospin effect [12J. For the contribution from target fragmentation this ratio should 
be constant into the backward hemisphere. For the determination of the feed-down correction 
the corresponding yields have to be folded with the on-vertex baryon reconstruction efficiency 
which reaches large values in the far backward hemisphere. The resulting correction in percent 
of the total proton yield drops however quickly below xp < -0.2 due to the decrease of the A 
cross section relative to protons and due to the fact that the important baryon contribution from 
intranuclear cascading has no hyperon content. The numerical values in percent of the baryon 
yields are shown in Fig. |23]as a function of x^. 




Figure 23: Feed-down correction as a function of x^? at different px values for a) protons, b) 
anti-protons and c) neutrons, in the latter case integrated over 

In comparison to the pion data Pl this correction reaches considerable values of up to 
20% for the anti-protons and therefore constitutes, together with the absolute normalization, the 
most important source of systematic uncertainty. 

5.4 Re-interaction in the target 

The carbon target has an interaction length of 1.5%, which corresponds to about half of 
the length of the hydrogen target used in p-i-p collisions. The expected re-interaction correction 
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is therefore smaller than +0.5% in the forward and -2% in the backward hemisphere. The values 
obtained in Q have therefore been downscaled accordingly. 

5.5 Absorption in the detector material 

The absorption losses in the detector material are equal to the ones obtained in [[3]|. 
Baryons in the newly exploited region in the far backward direction feature short track lengths 
in the first NA49 TPC detector only and are not affected by any support structures; hence the 
corresponding corrections are below the 1% level as shown in Fig. [241 




-0.5 0.5 1 -0.2 



Xp Xp 

Figure 24: Detector absorption correction as a function of at different px values for a) protons 
and b) anti-protons 



5.6 Binning correction 

The correction for finite bin width follows the scheme developed in [IJ using the lo- 
cal second derivative of the particle density distribution. This correction stays, despite of the 
rather sizeable bins used in some areas of the p-i-C data, generally below the ±2% level, being 
negligible in the xp co-ordinate for protons due to their rather flat dn/dxp distribution. Two 
typical distributions of the px correction for protons and for the xp correction for anti-protons 
are shown in Fig. [25] both for the nominal bin width of 0.1 GeV/c in px and 0.05 in xj? and for 
the actually used bin widths. 




[GeV/c] Xp 

Figure 25: Binning correction a) in pt for protons dX xp = Q and h)'m xp for anti-protons at 
Pt = 0.3 GeV/c. The values for the nominal bin widths of 0.1 GeV/c in px and 0.05 in xp are 
indicated by crosses, for the actually used bin widths by open circles 
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5.7 Systematic errors 



An estimation of the systematic errors induced by the absolute normalization and the 
applied corrections is given in Table \T\ For the proton and anti-proton extraction in the newly 
exploited backward regions an additional systematic error due to particle identification is indi- 
cated. An upper limit of 7.0% (8.5%) for protons (anti-protons) results from the linear addition 
of the error sources, increasing to 10% (14.5%) in the backward region. Quadratic summation 
results in the corresponding values of 3.7% (4.2%) and 4.7% (7.3%). The distribution of the 
numerical values of the corrections in all measured bins is shown in Fig.[26]for protons. 
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Table 1 : Systematic errors 
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Figure 26: Distribution of proton corrections for a) target re-interaction, b) empty target, c) 
trigger bias, d) absorption, e) feed-down, f) binning, g) total 



6 Results on double-differential cross sections for p and p 
6.1 Data tables 

The binning scheme presented in Sect. [3] results in 491 and 121 data values for protons 
and anti -protons, respectively. These are presented in Tables |2] and |3l 
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Table 2: Invariant cross section, f{xF,PT), in mb/(GeV^/c^) for protons in p+C collisions at 
158 GeV/c beam momentum. The relative statistical errors, A/, are given in %. The systematic 
errors are given in Table [T] 
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Table 3: Invariant cross section, f{xF,PT), in mb/(GeV^/c^) for anti-protons inp-i-C collisions at 
158 GeV/c beam momentum. The relative statistical errors, A/, are given in %. The systematic 
errors are given in Table [T] 

6.2 Data interpolation 

As in the preceding publications concerning p+p and p+C |^ interactions, a two- 
dimensional interpolation is applied to the data which reduces the local fluctuations given by 
the errors of the data points by a factor of 3-4. As there is no possibility to describe the de- 
tailed xp and pt distributions by simple functions and as any algebraic approximation tends 
to dilute the data quality by introducing systematic biases, the chosen interpolation relies on a 
multi-step recursive method using eyeball fits. The quality of the procedure may be controlled 
by plotting the differences between data points and interpolation, normalized to the statistical 
errors. The resulting distribution should be a Gaussian centred at zero with variance unity. This 
is demonstrated in Fig.|27]for protons and anti-protons. 




A/A f A/A f 



Figure 27: Plots of the normalized difference between data and interpolation for a) protons and 
b) anti-protons 



20 



6.3 Dependence of the invariant cross sections on pt and xp 

The distribution of the invariant cross section as a function of px is shown in Fig. [28] 
for protons and anti-protons at negative xp and in Fig. |29]at positive xp, indicating the data 
interpolation by full lines. For better visibility successive values mxp have been multiplied by 
a factor 0.5. 




0.5 1 1.5 2 0.5 1 1.5 2 

[GeV/c] [GeV/c] 

Figure 28: Invariant cross sections and data interpolation (full lines) as a function of pt at 
fixed Xi? < for a) protons and b) anti-protons produced in p-i-C collisions at 158 GeV/c. The 
displayed cross sections have been multiplied by a factor of 0.5 for successive values ofxp for 
better visibility 
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Corresponding xp distributions are presented in Fig. [30] for protons and in Fig. 1311 for 
anti-protons. 

The shape of the px distributions resembles, for x^? > (Fig. [291), the one measured 
in p+p interactions [3] including details of the deviation from either exponential or Gaussian 
shape. In backward direction. Fig. [28l a more complex behaviour develops with a steepening 




0.5 1 1.5 2 0.5 1 1.5 2 

[GeV/c] [GeV/c] 

Figure 29: Invariant cross sections and data interpolation (full lines) as a function of pt at 
fixed Xi? > for a) protons and b) anti-protons produced in p-i-C collisions at 158 GeV/c. The 
displayed cross sections have been multiplied by a factor of 0.5 for successive values of xp for 
better visibility 
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up at Xi? < -0.4. The basic asymmetry of the p+C interaction is more directly visible in the xp 
distributions of Figs. [30] and 31. If at highp^ factors of 1.6-2 are typical between the backward 
and forward proton yields at xp = 0.5-0.7, these factors grow to about 3-5 at pp = 0.1 GeV/c. 
This allows a first view at the composition of p+A collisions from projectile fragmentation 




-1 -0.5 0.5 1 

Figure 30: Invariant cross sections as a function of Xi? at fixed pp for protons produced in p-i-C 
collisions at 158 GeV/c. The px values are to be correlated to the respective distributions in 
decreasing order of cross section 
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in forward direction and target fragmentation as well as intranuclear cascading in backward 
direction. A quantification and separation of these three basic ingredients will be performed in 
Sects. [I2l[l3] and [14] below. For anti-protons the xp distributions clearly peak at negative xp- 
The forward-backward asymmetry is, at \xf\ = 0.2, to first order p-r independent and on the 
order of 1.6-1.9 which is well above the asymmetry of protons at this xp indicating important 
effects from isospin and baryon number transfer. 



CO 

O 

> 
CD 

CD 



10 



10" 



10" 



10 



-3 



"1 I r 



1 \ r 



pC^pX 




-0.2 







0.2 



X, 



Figure 31: Invariant cross sections as a function of xp at fixed pr for anti-protons produced in 
p-i-C collisions at 158 GeV/c. The px values are to be correlated to the respective distributions 
in decreasing order of cross section 
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6.4 Rapidity and transverse mass distributions 

The rapidity distribution for protons at fixed px presented in Fig. [32l extends to the 
kinematic limit in forward direction and to -2.6 units in the target hemisphere. Again a clear 




y 

Figure 32: Invariant cross sections as a function of y at fixed pr for protons produced in p+C 
collisions at 158 GeV/c. The pr values are to be correlated to the respective distributions in 
decreasing order of cross section 
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view of the asymmetry of the p+C interactions increasing with decreasing pt is evident. 

The rapidity range of anti-protons is limited by statistics to -1.4 to -1-1.8 units. The cor- 
responding distribution as a function of y for fixed is shown in Fig. [331 




Figure 33: Invariant cross sections as a function of y at fixed px for anti-protons produced in 
p-i-C collisions at 158 GeV/c. The px values are to be correlated to the respective distributions 
in decreasing order of cross section 
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Transverse mass distributions at ?/ = 0, with uit 



+ p^, are presented in Fig. 



together with the local inverse slopes. A situation very similar to p+p collisions emerges with 
a non-exponential behaviour and inverse slope parameters varying strongly with the transverse 
mass. 
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Figure 34: Invariant cross section as a function of uit — mp for a) protons and b) anti-protons. 
Panels c) and d) give the inverse slope parameters of the itlt distributions as a function of 
niT — mp for protons and anti-protons, respectively. The full lines represent the results of the 
data interpolation 



6.5 Comparison to other experiments 

As stated in Sect. [21 there are only two experiments providing double differential data in 
the SPS energy range. The data set [7] which is disjoint from the NA49 phase space coverage 
will be discussed in detail in the next Sect. |7J As far as the Fermilab data of Barton et al. [9l 
are concerned there are 10 data points for protons and 4 data points for anti-protons available in 
overlapping phase space ranges. The situation for protons is shown in Fig.[35]where the invariant 
cross sections [9] at px = 0.3 and 0.5 GeV/c are shown together with the data interpolation of 
the NA49 experiment as a function of xp. 

The comparison between the two experiments reveals very sizeable systematic offsets 
with an average of -1-21% or +6 standard deviations averaged over all data points. This complies 
with the comparison of pion yields JH with an average of +25% or -1-3.6 standard deviations. 
The good agreement between the particle ratios (K+-i-7r+)/p demonstrated in Sect. l4.1.4l over the 
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Figure 35: Invariant proton cross sections from for pj- = 0.3 and 0.5 GeV/c together with the 
data interpolation of NA49 (full lines) and the interpolation multiplied by 1 .21 (dashed lines) 



full comparable range oixp and "Pt speaks indeed for a normalization problem as the origin of 
the discrepancies which are not visible in the results of p+p interactions Q. 

For anti-protons the situation is considerably less clear due to the very large statistical 
errors of the data [9| indicating only upper limits for some of the measurements. Fig. [361 gives 
the xp dependence of these data for = 0-3 GeV/c together with the NA49 interpolation (full 
line) which has been partially extrapolated using an exponential shape. The large upward shift 
of the Fermilab data of about a factor of 2 (dashed line in Fig. 136]) is compounded by the fact 
that due to the lower beam momentum of 100 GeV/c these data should be expected to be about 
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Figure 36: Invariant anti-proton cross sections from |l9l for = 0.3 GeV/c as a function oixp 
together with the data interpolation of NA49 (full line) which has been extrapolated Xa xf = 0.6 
with an exponential function (dotted line). The NA49 reference multiplied by a factor of 2 is 
indicated by the dashed line 
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30-40% below the NA49 cross sections 

In this context it should be mentioned that already in p+p collisions the Fermilab data 
inn were high by about +25% for anti-protons. 

7 Data extension into the far backward direction 

As the backward acceptance of the NA49 detector is limited to the ranges of xf > -0.8 
to xp > -0.5 at low and high p^, respectively, it is desirable to extend this coverage into the 
far backward region to xp values down to and below -1. A detailed survey of existing experi- 
ments in the backward direction of p-i-C interactions has therefore been undertaken and is being 
published in an accompanying paper [6J. This survey establishes the detailed s-dependence of 
measured cross sections for beam momenta between 1 and 400 GeV/c, for lab angles between 
10 and 1 80 degrees, and for lab momenta between 0.2 and 1 .2 GeV/c. It shows in particular that 
the baryonic cross sections in the SPS energy range from about 100 to 400 GeV/c beam mo- 
mentum may be regarded as s-independent within tight systematic limits of less than a couple 
of percent. This allows the combination of the extensive data set of the Fermilab experiment |I3 
at 400 GeV/c beam momentum and lab angles between 70 and 160 degrees with the NA49 data 
which span the angular range up to 40 degrees. 

The relevant kinematic situation is presented in Fig. [37] where lines of constant piab and 
6iab are shown in the Xp/pr plane. 




-3 -2 -1 

Figure 37: Kinematics of fixed piab and Giab in the xplpp plane 

In the necessary transformation between the lab and cms frames involved in Fig. [37] 
there is very little difference between the beam momenta of 158 and 400 GeV/c. This is shown 
in Fig. [38] which gives the difference in a; as a function of piab for the angular range between 
70 and 160 degrees. 

One example of the apparent s-independence of the backward baryon yields is shown in 
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Figure 38: Differences inxp: Axp resulting from the transformation into the cms 

system between beam momenta of 158 and 400 GeV/c as a function of piab- The differences are 
independent of ©lab in the range 70 < ©lab < 160 degrees 



Fig. [39] which compares data at 6iab = 162 and 160 degrees for 8.5 [fT3ll and 400 GeV/c beam 
momentum Q. 




Figure 39: Invariant proton cross sections lfT3l at 0iab = 162 degrees and 8.5 GeV/c beam mo- 
mentum as a function of piab in comparison with the data from [7 J at 160 degrees and 400 GeV/c 
beam momentum. The full line represents the interpolation of the Fermilab data 

Due to the flat angular distribution at Oiab around 160 degrees, see Fig. |4TJ the small 
angular difference in angle between the two measurements has negligible influence on the cross 
sections comparison. Also the SPS measurements of grey protons by Braune et al. [[T4ll show 
no dependence on beam momentum in the range from 50 to 150 GeV/c over the complete 6iab 
range from 10 to 159 degrees. 

7.1 NA49 results at fixed ©lab and piab combined with the data from El 

The kinematic situation presented in Fig.[37]shows that the NA49 acceptance allows the 
measurement of proton yields as a function of pi^h up to 9iab = 40 degrees. The corresponding 
data values are tabulated in Table SI 
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Table 4: Invariant proton cross sections measured by NA49 at fixed values of ©lab between 10 
and 40 degrees as a for pi^h values between from 0.3 to 1 .6 GeV/c. The relative statistical errors, 
A/, are given in % 

These data are shown in Fig. I40l together with the Fermilab data [|7|| as a function of piab. 
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p,^JGeV/cl 

Figure 40: Invariant proton cross sections from NA49 and |I3 at lab angles between 10 and 160 
degrees as a function of piab- The data interpolation at fixed angle is given by the full lines 

Evidently the two data sets are complementary and offer for the first time an almost 
complete angular coverage of the backward proton production with double differential cross 
sections. 

A two-dimensional data interpolation has been performed as shown by the full lines in 
Fig. SOI This allows to produce the combined angular distribution as a function of cos(6iab) for 
fixed values of pub presented in Fig. Hi] 
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COS i&J 

Figure 41 : Invariant proton cross sections of the combined NA49 and Fermilab data as a function 
of cos(0iab) for fixed values of piab between 0.2 and 1.4 GeV/c. The data interpolation is shown 
as full lines 



The proton density distributions dup/dpiah derived from the interpolated invariant cross 
sections are shown as a function of piab in Fig. |42l surface normalized. 



Q. 




Figure 42: Surface normalized proton density distributions dnp/dpi^b as a function of piab for 
fixed Giab- The bubble chamber data [15] are shown as full circles 



These distributions are closely similar for 160 > 6iab > 1 1 8 degrees and develop a tail to 
large piab values for angles smaller than 90 degrees indicating increasing contributions from the 
fragmentation of the participant nucleons hit by the projectile. The "grey" proton momentum 
distribution measured in the EHS rapid cycling bubble chamber [15] is positioned close to the 
©lab = 90 degrees result. In this case a strong momentum cut is introduced by requesting bubble 
densities at 1.3 minimum ionizing rejecting most of the faster forward region. 
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The drip/dpiah distributions shown in Fig.l42lmay be integrated over piab resulting in the 
proton densities dn/dVt shown in Fig. |43]as a function of cos (6iab). Here the integration has 
been limited to piab < 1.6 GeV/c as at low angles the target fragmentation contribution will 
create a divergent behaviour and since the comparison data do not contain this component. 

c 

T3 



10-^ 



10-2 

-1 -0.5 0.5 1 
COS(0|ab) 

Figure 43: Proton density dn/dVt as a function of cos (©lab). Full circles: integrated combined 
data from NA49 and [7]. Open circles: measurement by Braune et al. [14J 

The direct measurements from fT4l| also presented in Fig. |43] show systematic deviations 
both at forward and at backward angles. As this experiment uses an energy loss measurement 
with variable threshold it is not clear to which extent it represents identified proton yields. 
The authors in fact explain that a contribution from "evaporation" particles ("black tracks" in 
emulsion work) cannot be excluded. Such a contribution would typically be characterized by a 
flatter angular distribution as compared to protons, see also the discussion on light ions in Sect. [8] 
showing very sizeable d/p ratios. As the d/p and t/p ratios decrease steeply with increasing xp, 
see Fig.|45l the yield measured by lfT4l would increase at ©lab > 90° and decrease towards small 
angles with respect to the one of identified protons. In addition, a dE/dx cut-off in the detector 
of [11411 will reduce the contribution of fast protons from target fragmentation at small angles. 

7.2 Double differential cross section f (xp, pt) as a function of xp 

The invariant cross sections measured in piab and 0iab may be transformed into the x p/pr 
variables following the kinematics shown in Fig.[37l The corresponding distributions at fixed 
Pt are presented in Fig. l44l which shows the xp range down to -2.0. Both the NA49 data (full 
circles) and the Fermilab results (open circles) are plotted. The full lines in Fig. l44l represent 
the interpolation of the two data sets also covering the non-measured angular region between 
40 and 70 degrees, see Fig. |4T] The thin lines indicate the position of the lab angles of 10, 30, 
and 50 degrees as well as the five angles measured by [7]. 

Fig. l44lrepresents one of the main results of this paper. It shows for the first time a com- 
plete coverage of the baryonic phase space in p-i-A collisions, from a;^ = -2 up to the kinematic 
limit for the projectile fragmentation atxp ~ -i-l. Several features merit comment: 

- the invariant proton cross sections extend far below the kinematic limit for target frag- 
mentation at Xi? = -1 




33 




-2-10 1 

Figure 44: Invariant cross sections at fixed as a function oixp- Full circles: NA49 data, open 
circles: data from [TJ. The thin lines show the cross section at fixed angles of 10", 30° and 50° 



- there is no indication of a diffractive structure with a peak xp ~ - 1 as it would be 
expected from the prompt fragmentation of the hit target nucleons - on the other hand, 
diX. xp > +0.9 there is a diffractive peak in the projectile fragmentation region, see |l3l 
for comparison with p-i-p interactions 

- the backward cross sections peak at x^? ~ -0.9, not dX xp ~ -1 indicating a sizeable 
longitudinal momentum transfer in the nuclear fragmentation region 

- at low pt or low transverse momentum transfer however, the lines of constant lab angle 
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are compatible as expected with a convergence towards a;^ = -1 

A detailed discussion of these features including the de-composition of the measured 
proton yields into the basic components of projectile, target and nuclear fragmentation, is pre- 
sented in Sects. [T3] and [141 

8 Light ions: deuterons and tritons 

As shown in Sect. HI Fig. [51 the particle identification via energy loss measurement in 
the NA49 TPC system also allows the extraction of deuteron and triton yields. The accessible 
kinematic region covers lab momenta from the detector acceptance limit at about 0.25 GeV/c 
up to the crossing of the energy loss distributions with the ones for electrons. Fig. [3l at about 
Piab = 2 GeV/c for deuterons and 3 GeV/c for tritons. 

As for the proton cross sections, this range is complementary to the Fermilab experiment 
[[H which gives light ion cross sections in the lab angular range from 70 to 160 degrees, at pi^^ 
from 0.5 to 1.3 GeV/c for deuterons and from 0.7 to 1 .3 GeV/c for tritons. The NA49 data offer 
the advantage of reaching low pt down to 0. 1 GeV/c and of covering the forward region from 
©lab = 40 degrees down to about 3 degrees. This allows for the first time to trace the extension 
of nuclear fragmentation into light ions towards the central region of particle production. 

8.1 Ion to proton ratios 

In order to clearly bring out the last aspect mentioned above the deuteron and triton 
yields are given here as ratios to the proton yields in each bin of piab and 6iab, 
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These density ratios are given as functions of xp and pt using the proton mass in the 
transformation from lab to cms variables. They thus give directly the relative contribution of the 
light ions with respect to protons in the xplpT bins shown in Fig. [2| and Table [2l In forming the 
ion/proton density ratios, most of the data corrections. Sect. [5l drop out with the exception of 
absorption which is small but increased for ions, and of hyperon feed-down which is of course 
only applicable to protons. The resulting ratios are shown as functions oi xp and px in Table [51 
for deuterons and in Table [6l for tritons. 
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Table 5: d/p density ratios Rd{xF,PT) and their statistical errors ARd in percent as a function 
of xf and pr using proton mass in the transformation from lab to cms system 
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Rt{xp,pT], ARt 


Pt\xf 


-0.8 


-0.75 


-0.7 


-0.65 


-0.6 


-0.55 


-0.5 


-0.45 


-0.4 
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-0.3 


0.1 


1.22 


0.1 




1.74 0.9 


0.79 0.5 




1.18 0.7 


1.00 0.3 


0.44 0.2 


0.31 0.2 


0.45 0.3 
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1.12 


0.4 
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0.34 0.3 
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0.44 0.2 
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0.15 0.1 


0.28 0.2 


0.41 0.3 
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0.82 0.4 
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0.69 0.3 


0.99 0.5 


0.43 0.3 


0.73 0.4 


0.40 0.3 


0.47 0.2 


0.5 










2.08 0.8 


1.16 0.5 


0.33 0.3 


0.67 0.4 


0.38 0.3 




0.20 0.2 


0.36 0.2 
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0.45 0.4 




0.48 0.4 









Table 6: t/p density ratios Rt^xp, Pt) and their statistical errors ARt in percent as a function of 
xp and Pt using proton mass in the transformation from lab to cms system 



In these data tables the cut-offs imposed by the NA49 acceptance at 6iab ~ 40 degrees 
below xp ^ -0.5 and by the upper limit on pi^h imposed by the energy loss measurement at 
> -0.5 are discernible. In addition the fast decrease of the ratios towards higher xp limits 
the extraction of tritons due to the low overall event statistics. 



8.2 Comparison to the Fermilab data llsll 

The density ratios Rd{xp,pT) and Rt{xp,pT) are presented in Fig. |45]as a function of 
xf in comparison to the Fermilab data which are available above pp ~ 0.3 GeV/c. 
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Figure 45: Deuteron and triton to proton density ratios Rd and Rt as a function of xp for fixed 
values of pp between 0.1 and 0.9 GeV/c. The full lines give the two-dimensional data interpo- 
lation established for Rd, the broken lines are the same multiplied by the suppression factors 
shown in Fig. |46] 
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As for the proton data, the good consistency with the Fermilab data in the pr range 
above 0.3 GeV/c is confirmed for Rd by the two-dimensional data interpolation shown as the 
full lines in Fig. |45l Rt shows a similar consistency with [8] over the smaller available region 
of comparison due to the higher piab cutoff in the Fermilab data. With respect to R^,, Rt is 
suppressed by a factor of 0.1 to 0.2 depending on px but within statistics independent of xp- 
This is demonstrated by the broken lines in Fig. |45] which represent the interpolation of Rd with 
the suppression factors Rt/Rd^^'^"^ given as a function of pr in Fig. |46l 

o 0.2 1 

e- 
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0.2 0.4 0.6 0.8 1 
P, [GeV/c] 

Figure 46: Triton suppression factor RtlR^f^^"^ as a function of pt 



8.3 xf and px dependences 

The xp dependences of the density ratios presented in Fig. |45] show a rather complex 
behaviour. Evidently deuterons and tritons, if seen as nuclear fragments as opposed to coalesc- 
ing from produced baryons, reach far towards central production. In consequence the separation 
of nuclear fragmentation and coalescence is not an easy affair especially as fragmentation prod- 
ucts, as shown by the tentative extrapolation of the data interpolation towards x^r = -0.1, might 
well represent a contribution in the sub-percent range even dXxp = and already for the light 
Carbon nuclei at SPS energy. The pr dependence of Rd at fixed is as well non-trivial, as 
shown in Fig.l47l 




[GeV/c] 

Figure 47: Rd as a function of pt for fixed values of xp between -1.4 and -0.2 

Evidently Rd tends to decrease with increasing pp for > -0.8, this tendency inverting 
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itself for < -0.8. The RtlRa ratio on the other hand clearly increases with pr at all available 
xf values, see Fig.l46l 



9 Particle ratios 

The already published data on p+C [|4l[5l and p+p [[T]43l interactions allow for a very 
detailed study of particle ratios. In a first overview baryonic ratios will be investigated in this 
section, both concerning anti-proton/proton ratios in elementary and nuclear collisions. 



m = f^'i^F,PT)/f^n^F,PT) 

and baryon density ratios directly comparing p+C and p+p reactions, 



(6) 
(7) 



R, = ffixF,PT)/f^n^F,PT)- 



(8) 
(9) 



Fig. |48] shows R^^ (closed circles) and R^^ (open circles) together with the correspond- 
ing ratios of the data interpolation (full and broken lines) as a function of for several bins of 
Pt- 
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Figure 48: Rpp as a function of for 8 values of pt comparing p+C (closed circles) and p+p 
(open circles) interactions. The corresponding ratios of the data interpolations are shown as full 
lines (p+C) and broken lines (p+p) 



This comparison reveals that the ratios are within errors equal in the projectile hemi- 
sphere xp > 0.1 with the exception of the two highest pr bins. This means that the transfer of 
projectile baryon number ("stopping") towards the central region is equal for protons and for 
anti-protons. In the backward hemisphere there are distinct differences between p+C and p+p 
interactions at low to medium px- This is a result of isospin effects between the isoscalar C 
and the p targets. In fact it is known that the anti-proton yields increase and the proton yields 
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decrease in neutron fragmentation lfT2]| with both effects increasing R^^ compared to R^. The 
clean extraction ofp and p yields in the far backward direction, a; < -0.4, see Sect.|4]and Fig.[5l 
allows the extension of R^^ to the xp range -0.7 < xp < -0.5. This is presented in Fig. 
the pt interval OA < pt < 0.4 GeV/c. 
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Figure 49: R^^ in the intervals -0.7 < xp < -0.5 and 0.1 < pt < 0.4 GeV/c together with R^^ 
(full line), (broken line) m and IJR^l (dotted line) 



The broken line represents R^ [3J, the full line the measured R^^ and the dotted line R^ 
multiplied by a factor 1.7 which stems from the expected increase in p yield (factor 1.33, [,12J) 
and the decrease in p yield (factor 1.3, O) from the isoscalar C nucleus. The fact that the 
measured backward ratio follows closely the expectation from target fragmentation shows that 
there is no p production from nuclear cascading. This is in agreement with the upper limit of 
^ given in [16] for p-i-^^^Ta at 90 degrees laboratory angle. 
In contrast to the p/p ratios discussed above, the baryon yields proper exhibit an impor- 
tant evolution when passing from p-i-p to p-i-C interactions. This is presented as functions of xp 
for fixed values of p^ in Fig. [50]for Rp and in Fig. [5T]for _Rp. 
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Figure 50: Rp as a function of xp for fixed px values between 0.1 and 1.3 GeV/c. The full lines 
give the ratios for the corresponding data interpolations 
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Figure 51: i?p as a function of for fixed pr values between 0.1 and 1.3 GeV/c. The full lines 
give the ratios for the corresponding data interpolations 



Looking first at Rp, Fig.[50l where the data give access to a very wide range ofxp, three 
distinct zones may be distinguished. 

A first area is governed by the projectile fragmentation, xp ^0.1. Here the transfer 
of baryon number from the forward to the central region is clearly evident, with a suppression 
of the far forward yield by a factor of 0.6 at low to medium pt followed by a steady increase 
of proton density with increasing px until the densities exceed the p+p values at all xp for 
Pt > 1 GeV/c. 

A second area is characterized by the fragmentation of the target nucleons hit by the 
projectile, -0.5 < xp < 0. Here the pile-up of produced hadrons in the central area due to 
the 1.6 mean projectile collisions in the C nucleus is visible, followed by a steady decrease of 
proton density in the backward direction due to the contribution of neutron fragmentation. 

A third component is characterized by the steep increase of proton density for xp < - 
0.2, at low to medium pp. This region is governed by the intranuclear cascading of baryons 
("grey protons") in the lab momentum range up to about 2 GeV/c. This component dies out as 
expected with increasing transverse momentum. 

Due to the smaller xp range experimentally available for the anti-protons, the corre- 
sponding ratio Rp is more restricted. It has however been demonstrated (Fig. |49l) that here no 
cascading component exists. As argued above, the evolution of the projectile fragmentation 
(xp ^0.1) follows closely the one for protons including the increase with transverse momen- 
tum. As far as the target fragmentation is concerned, the increase of anti-proton density beyond 
the pile-up from the 1.6 average projectile collisions to about a factor of 2 is expected from 
isospin symmetry. 

A further, more detailed argumentation concerning this phenomenology can be found in 
the discussion Sect. [T2l below. 

10 Integrated data 

10.1 Pt integrated distributions 

The Pt integrated non-invariant and invariant yields are defined by: 



40 



dn/dxp = vr/cTinei ■ -\/s/2 ■ j f / E ■ dp^ 

f ■ dvl (10) 



dn/ dy = vr/ cTinei ' j f ' ^Pr 



with f = E ■ d^a/dp^, the invariant double differential cross section. The integrations are 
performed numerically using the two-dimensional data interpolation (Sect. [6^ . Tables |7] and [8] 
give the numerical values and the first and second moments of the px distributions as functions 
of xf and rapidity for protons and anti-protons, respectively. The relative errors quoted are 
given in percent using the full statistical errors of the measured points. They thus present upper 
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-0.45 


32.8 


1.8 
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Table 7: p^ integrated invariant cross section F [mb-c], density distribution dn/dxp, mean 
transverse momentum {pj) [GeV/c], mean transverse momentum squared (p^) [(GeV/c)^] as 
a function of Xi?, as well as density distribution dn/dy as a function of y for p. The statistical 
uncertainty A for each quantity is given in % as an upper limit considering the full statistical 
error of each measured pr/xp bin 
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limits as the data interpolation inxp will reduce the statistical fluctuations for the pr integrated 
quantities. 
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Table 8: pr integrated invariant cross section F [mb-c], density distribution dn/dxp, mean 
transverse momentum {px) [GeV/c], mean transverse momentum squared (p^) [(GeV/c)^] as 
a function of xp, as well as density distribution dn/dy as a function of y for p. The statistical 
uncertainty A for each quantity is given in % as an upper limit considering the full statistical 
error of each measured pt/x f bin 

Concerning the extension of the kinematic coverage into the far backward direction 
for protons, xp < -0.8, see Sect. |7l additional integrated quantities are given in Table |9] for 
-2.0 <xf< -0.9. Again, the data interpolation is used to obtain the numerical values. 
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0.400 


2.0 


0.211 


3.0 


-1.1 


71.0 


3.0 


0.890 


3.0 


0.384 


2.0 


0.197 


3.0 


-1.0 


88.8 


3.0 


1.224 


3.0 


0.374 


2.0 


0.189 


3.0 


-0.9 


93.0 


3.0 


1.422 


3.0 


0.378 


2.0 


0.194 


3.0 



Table 9: px integrated invariant cross section F [mb-c], density distribution dn/dxp, mean 
transverse momentum {px) [GeV/c], mean transverse momentum squared (p^) [(GeV/c)^] as a 
function of Xi;' for protons in the far backward region. The estimated relative statistical uncer- 
tainties A are given in % 

In this region the estimation of the statistical uncertainties is non-trivial as the measured 
data have been transformed from the laboratory co-ordinates (piab,0iab) to the cms quantities 
(xf,Pt) using a two-dimensional interpolation. The given errors have been obtained from the 
total number of measured protons at each lab angle [[H which varies from about 1 100 to 3600. 
The values given in Table |9] are therefore to be regarded as upper limits. 

The corresponding distributions for dn/dxF, F and dn/dy are shown in Fig. [52] for 
protons and in Fig.|53]for anti-protons. 

The pt integrated p/p ratios as a function oi xf are shown in Fig. [541 The first and 
second moments {pr) and {p^) are presented as a function of in Fig. [55] for protons and in 
Fig. [56] for anti-protons. 
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Figure 52: pt integrated distributions a) dn/dxF,h) F and c) dn/dy for protons as a function 
of xp and y, respectively 
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Figure 53: integrated distributions a) dn/dxF,h) F and c) dn/dy for anti-protons as a func- 
tion of xp and y, respectively 
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Figure 54: The pt integrated p/p ratios as a function of xp 
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Figure 56: a) mean px and b) mean p| for anti-protons as a function of xp 



10.2 Centrality dependence 

The detection of "grey" protons in the centrality detector of NA49 (W\ allows the study 
of particle densities as a function of the number n„^ey of protons in the lab momentum range 
of approximately 0.15 to 1.2 GeV/c. The number distribution dN/dn^^^y has been shown in 
BU to be a steep function of rigrey with only 20% and 5% of all events at ngrey = 1 and 2, 
respectively. The determination of double differential cross sections is therefore not feasible in 
this experiment. As already shown in [4J for pions, the extraction of pr integrated yields as 
a function of x^? by fitting the dE / dx distributions over the complete range of pt is however 
feasible in a limited range oixp- This range is determined by the variation of total momentum 
with Pt and extends from xp = -0.2 to -1-0.65 for protons, the upper limit being imposed by 
the progressive loss of acceptance at low px- Two samples with rigrey > 1 and rigrey > 2 have 
been selected. The resulting proton density distributions dn/dxp are given after correction in 
Table [TOl with the additional use of the complete data sample (see also [|4J) which allows for a 
precise control of eventual systematic effects. 

As shown in Fig. [571 the results for the total data sample are compatible within errors 
with the integration of the data interpolation, Tableland Fig.[52l 
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{dn/dxpf" 


A 




A 


[dn/dxpT"'-^ 


A 


-0.2 


0.797 


3.8 


1.210 


5.3 


1.388 


8.2 


-0.15 


0.858 


2.4 


1.225 


5.4 


1.461 


8.8 


-0.1 


0.827 


3.3 


1.230 


5.2 


1.401 


10.1 


-0.05 


0.852 


3.2 


1.247 


4.6 


1.373 


8.9 


0.0 


0.923 


3.0 


1.167 


2.7 


1.285 


6.0 


0.05 


0.867 


1.5 


1.143 


2.4 


1.201 


4.6 


0.1 


0.774 


1.4 


0.975 


2.2 


1.055 


4.0 


0.15 


0.725 


1.4 


0.824 


2.4 


0.903 


4.7 


0.2 


0.678 


1.5 


0.800 


2.4 


0.838 


4.8 


0.25 


0.682 


1.5 


0.758 


2.5 


0.738 


5.0 


0.3 


0.636 


1.6 


0.649 


2.7 


0.615 


5.5 


0.35 


0.594 


1.5 


0.597 


2.9 


0.541 


6.1 


0.4 


0.579 


1.6 


0.532 


3.0 


0.519 


6.2 


0.45 


0.510 


1.6 


0.462 


3.2 


0.426 


6.8 


0.5 


0.478 


1.7 


0.451 


3.3 


0.424 


7.2 


0.55 


0.441 


1.8 


0.384 


3.7 


0.360 


7.5 


0.6 


0.423 


2.8 


0.388 


5.7 


0.261 


13.7 


0.65 


0.383 


3.4 


0.312 


6.8 


0.298 


13.4 



Table 10: dn/dxp distributions for protons as a function of for the complete data sample and 
the selections Ugrey > 1 and ngj-ey > 2. The relative statistical errors are given in % 




Figure 57: dn/dxp for protons as a function of xp for the complete data sample and for the 
conditions a) rigrey > 1 and b) rigrey > 2. The broken line gives the density in p-i-p collisions |l3l 
and the full line the integration of the interpolated p-i-C data 

As seen from Fig. [57l a systematic and smooth variation of the proton density distribu- 
tions is evident when passing from p-i-p to minimum bias p-i-C and to centrality enhanced p-i-C 
interactions. A well defined cross-over with equal number density is visible dXxF = 0.25 ^ 0.3. 
Above this value the densities decrease progressively with increasing number of projectile col- 
lisions by up to a factor of 0.5 at the experimentally accessible limit of = 0.6. Below xf = 0.3 
this trend is inverted with a density increase of up to factors of 2.5 at negative xp- 
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A similar behaviour is also seen for anti-protons, although here only a reduced range in 
Xi? is accessible due to the limiting statistical uncertainties. This is shown in Fig. 158) where the 
relative difference in baryon density between the Ug^^y selected and minimum bias samples is 
given in % as a function of xp for protons (panel a) and anti-protons (panel b). 
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Figure 58: Relative difference A in % between ngrey selected and minimum bias samples as a 
function of xp for a) proton and b) anti-protons. The lines are shown to guide the eye. The lines 
in panel b) repeat those from panel a) 



The interplay between baryon number transfer from the projectile hemisphere and pile- 
up of baryons from the target fragmentation, including isospin effects, will be discussed in more 
detail in Sect. 12 of this paper. 

10.3 Neutron data 

The detection of forward neutrons in the Ring Calorimeter of NA49 [fTOll has been intro- 
duced and described in detail in tSJ. The neutron analysis concerning the separation of electro- 
magnetic and hadronic deposits, the veto against charged particles, the calorimeter calibration 
and the energy resolution unfolding may be directly applied to the p-i-C interactions. The cor- 
rections for acceptance and for feed-down from weak decays are shown in Figs. [2] and [231 
respectively. The empty target, trigger bias and re-interaction corrections are equal to the ones 
for protons |[3j. The contributions from Kl decay and anti-neutron production are obtained from 
the charged kaon data which are available from NA49 1,17,1 and from the isospin argumentation 
using the anti-proton data explained in [O. 

The resulting neutron densities dn/dxp are listed in Table [TT] together with their ratio to 
the results from p-i-p collisions lO. 

The pt integrated neutron density distribution dn/dxp as a function of a;^ is presented 
in Fig. [59k . Of particular interest in this context is the ratio between neutron densities in p-i-C 
and p-i-p interactions in comparison to the same ratio for protons as shown in Fig. [59b. 

Evidently protons and neutrons show within errors the same behaviour when passing 
from p-i-p to p+C collisions. The density ratio is equal to 1 at = 0.3, increasing to 1.3 at 
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[dn/dxpf^ /{dn/dxFf^ 


A 


0.1 


0.621 


21.4 


1.29 


5.0 


0.2 


0.482 


15.1 


1.18 


3.2 


0.3 


0.389 


13.6 


1.03 


3.1 


0.4 


0.303 


12.7 


0.93 


3.2 


0.5 


0.268 


12.8 


0.83 


3.6 


0.6 


0.221 


11.0 


0.75 


4.1 


0.75 


0.194 


11.7 


0.68 


5.2 


0.9 


0.128 


29.1 


0.59 


8.4 



Table 11: pt integrated neutron density distribution dn/dxp as a function of and the ratio 
of neutron densities in p+C and p+p interactions. The relative errors are given in %. They are 
governed by the systematic uncertainties quoted in Table [T] 
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Figure 59: a) pr integrated neutron density distribution dn/dxp and b) neutron density ratio 
between p+C and p+p interactions as a function oi xp. The full line shows the corresponding 
density ratio for protons 

xp = 0.1 and decreasing to 0.6 aXxp towards +1. The latter value corresponds to the expected 
fraction of single projectile collisions in p+C interactions derived from the nuclear density dis- 
tribution s. 

10.4 Comparison to other data 

Two data sets are available for comparison with the NA49 integrated data in the SPS 
energy range. The first set from the ACCMOR collaboration [18J gives the density distribution 
dn/dxp of protons for p+A collisions from Be to U nuclei at 120 GeV/c beam momentum. A 
second, very recent publication from the MIPP collaboration [fT9l provides neutron densities 
dn/dxp for p+C interactions again at 120 GeV/c beam momentum at the FERMILAB Main 
Injector. 

10.4.1 ACCMOR data ^ 

The ACCMOR collaboration [T8l measured proton densities dn/dxp with a proton 
beam at 120 GeV/c momentum at the CERN SPS for five different nuclei (Be, Cu, Ag, W, 
U) in the xp range from 0.07 to 0.6. These data may be interpolated from Be to C. This is 
demonstrated in Fig. [60] which gives the proton density as a function of mass number A for five 
xp values from 0.1 to 0.6. 
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Figure 60: Proton density dn / dxp as a function of mass number A for p+Be, p+Cu, p+Ag, p+W 
and p+U nuclei for different values oixp- The lines describe the data interpolation used for the 
determination of the corresponding p+C cross sections 

The data interpolation of dn/dxp to p+C interactions is shown in Fig.[6T]as a function 
of xp in comparison to the NA49 results, Table|7]and Fig. [52] above. 
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Figure 61: Proton density dn/dxp a function of Xi;' interpolated to p+C from the ACCMOR 
data. The full line through the ACCMOR results is drawn to guide the eye, the full line shows 
the NA49 result. The broken line gives the ACCMOR result after a tentative subtraction of 
proton feed-down from weak decays. The dotted line corresponds to the modification of the 
total inelastic cross section for p+Be collisions with 6%, see Fig. [62] 

The ACCMOR data show an upward deviation which increases from about atxp = 0.6 
to about 20% at the lower data limit of = 0.07. This deviation looks similar to the feed-down 
correction for protons from weak hyperon decays used for the NA49 data. Sect. 15.31 which 
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reaches 15% at small xp and decreases rapidly to at xp ~ 0.7. As such a correction is not 
mentioned in [[T81 and as the geometrical layout of the experiment is similar to the one of NA49 
as far as the position of the tracking elements is concerned (~ 3 m distance between target and 
first tracking station) it has been tentatively assumed that the on- vertex reconstruction efficiency 
of decay products might be similar in both experiments. Subtracting the feed-down correction 
from NA49 results in the broken line in Fig. [52l This line is about 8% above the NA49 results 
for Xi? < 0.4. At this point a look at the total inelastic cross sections as a function of A from 
different references lfT8ll20l . Fig. [621 shows that the cross section for p-i-Be given in [18] falls 
low by about 6% compared to the interpolation of other available data [[201 . A corresponding 
correction of the cross section is shown as dotted line in Fig. [61] 
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Figure 62: Data summary of measured total inelastic cross sections of p-i-A interactions as a 
function of A, [[T8[[20l 

In conclusion it may be stated that the ACCMOR results are compatible with the NA49 
results within the quoted systematic uncertainties. 

10.4.2 MIPP data 

Very recently new data on neutron densities dn/dxp have become available from the 
MIPP collaboration at the FERMILAB Main Injector for p-i-p and p-i-A interactions [19]. The 
MIPP neutron densities dn/dxp for p-i-C interactions at 120 GeV/c beam momentum are com- 
pared to the NA49 results in Fig. [63k . the densities for p-i-p collisions at 84 GeV/c beam mo- 
mentum in Fig. [63b . 

Important deviations are visible between the two experiments both for p-i-C and p+p 
collisions. This is quantified in Fig. [64[ where the ratio 

Rn = {dn/dxFr'''''/{dn/dxpf^^'' (11) 
between the two respective interpolations is shown as a function of xp- 
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Figure 63: Comparison of MIPP neutron densities dn/dxp to the NA49 results as a function 
of a;^, a) for p+C interactions and b) for p+p interactions. Interpolations through the respective 
data points are shown as full lines 




Figure 64: Density ratio Rn = (dn/dxp)^^^^ / (dn/dxp)^^'^^ as a function of xp, a) for p+C, b) 
for p+p interactions. The error bands shown are based on the statistical uncertainties of the data 
points 



The density ratios show a remarkable similarity between p+C and p+p interactions, vary- 
ing systematically from 1.4-1.6 at the lower xp limit to a maximum of 1.9-2.0 at a;^ ~ 0.4 
and decreasing to 0.5-1.0 at the upper kinematic limit. The total neutron yields in the forward 
hemisphere obtained by an extrapolation of the measured densities towards xp = as shown 
in Fig. [63] deviate by about 60% for both reactions, with the integrated MIPP neutron yield 
resulting in 1 neutron per event for p+p collisions at 84 GeV/c beam momentum. 

These large deviations pose a problem for baryon number conservation which may be 
verified for p+p interactions using published total yields for protons, neutrons, hyperons, pair 
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produced protons and pair produced hyperons. The available data are plotted in Fig. [65] as a 
function of y/s from ETll at 12 and 24 GeV/c beam momentum, It22ll23l at 32 GeV/c beam 
momentum, H24[|25ll at 69 GeV/c beam momentum and from NA49. 
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Figure 65: Total yields of protons, neutrons, hyperons, pair produced protons and anti-hyperons, 
[|2TI - l24]| . as well as the total number of net baryons (B-B) as a function of y/s between 4.9 and 
17.2 GeV. The integrated neutron yields from MIPP [[T9ll at 58 and 84 GeV/c beam momentum 
are indicated as open circles 

The large neutron yields published by MIPP are clearly out of proportion with respect 
to the available data and are incompatible with baryon number conservation. 

11 Two component mechanism of baryon and baryon pair production 

An important key to the understanding of hadronic interactions, in particular concern- 
ing the asymmetric p-i-A collisions, is provided by the fact that these may be separated into 
independent components connected to the fragmentation of the target and projectile systems, 
respectively. There are several experimental findings which prove this characteristics: 

- the absence, at SPS energy, of charge and flavour exchange 

- the absence of long-range two-particle correlations between the target and projectile 
hemispheres at \xf\ > 0.2 

- the independence of the target fragmentation on the nature of the projectile particle 
(hadronic factorization) 

There is however an overlap of target and projectile fragmentation in the central region at 
|a;ir| < 0.2, where this xp limit depends on the mass of the secondary particle considered and on 
the transverse momentum [5|. This has been experimentally verified by the measured forward- 
backward multiplicity correlations which extend - as they are governed by pion production - 
to l^i?! < 0.1. A detailed study of extent and shape of this overlap region has been performed 
for pions with the NA49 detector [5J using the isospin symmetry of the tt+Ztt^ ratio in 7r+-i-p 
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and 7r~+p interactions at the SPS. It has also been shown in ^ that the range and shape of the 
pionic overlap is equal for p+C and p+p collisions. In the following discussion the overlap from 
projectile to target hemisphere and vice versa will be called "feed-over". 

Concerning baryons the situation is somewhat more complicated as a clear distinction 
between pair-produced baryons and "net" baryons has to be made. Here the "net" baryon yield is 
defined as the difference between the total yield of a given baryon species and the yield of pair- 
produced baryons of the same species. Whereas the sum of all net baryon yields is constrained 
by baryon number conservation, the production of baryon/anti-baryon pairs shows a strong y/s 
dependence. The yield rises sharply from threshold at about ^/s = 6 GeV through the SPS 
energy range and flattens towards collider energies in a way characteristic of the production 
of heavy hadronic systems. The cross section is centred close to a;^ = , the central p/p ratio 
reaching about 23% at SPS energy in the minimum bias p-i-C interactions considered here, see 
Sect. \T0\ For ease of distinction, the pair produced baryons will be denoted in the following 
sections by: 



pair produced protons = p 
pair produced neutrons = n 

The situation is further complicated by the fact that baryon/anti-baryon pairs are not 
limited to the isospin combination but form an isospin triplet. In the case of protons and 
neutrons this results in the dependence on the isospin-3 component shown in Table [T2l 



-1 1 



baryon pairs 



pn 



PP 
nn 



np 



Table 12: Isospin structure of baryon pair production 



This phenomenology has been studied by NA49 using the comparison of p-i-p and n-i-p 
interactions. Here it was shown that the anti-proton yield increases (by a factor of about 1.6) 
when exchanging the / = -1-1/2 against an / = -1/2 projectile lfT2l . thereby enhancing the /s = -1 
with respect to the Is = -i-l combination. 

Both the baryon pair and the net baryon feed-over will be studied in the following sec- 
tions. In fact the extent and the shape of the corresponding overlap functions will be determined 
in a strictly model-independent way essentially based on measured quantities and isospin in- 
variance. In a first part, the symmetric p-i-p interactions will be exploited. In a second part it will 
be shown that the same overlap functions apply to p-i-C collisions. 



11.1 Two component mechanism in p+p interactions 

A large sample of 4.8 million inelastic events is available from NA49 [[I]-l3l with both 
proton/anti-proton and neutron identification in the final state. This allows for the selection of 
sub-samples of events with defined net baryon number either in the projectile or in the target 
hemisphere by tagging baryons at sufficiently large \xf\- Sufficiently large means in this context 
xp values where the yield of pair produced baryons is low enough to ensure negligible back- 
ground. For protons this condition is fulfilled for \xf\ > 0.35 where the p/p ratio is less than 
0.5% [l3l. The following ranges of xp have been used for net proton and neutron selection: 
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Projectile hemisphere: protons 

neutrons 

Target hemisphere: protons 



0.35 <xf < 0.5 
0.5 <xf < 0.7 
-0.75 <xf < -0.6 



(13) 
(14) 
(15) 



These ranges are given by the constraints of acceptance and proton identification via 
dE/dx, see Sect. HI and by the limited energy resolution of the hadron calorimetry for neutrons. 
In the following the measured baryonic double and single differential densities obtained with 
net baryon constraint, 



p\xf,pt) = -r—j— 

dXFdpT 

, (Id) 

dn 

Pint 



dxF 

will be described by their ratio to the corresponding inclusive yields, 



P'{xf,Pt) 

Pint(^F) 

by 

R^ixFyPr) 

The double differential cross sections are obtained, due to limited statistics, in the range 
0.1 < pt < 0.7 GeV/c, and the resulting ratios are averaged over this pr range resulting in the 
ratios: 



P%xf,Pt) 
P'{xf,Pt) 



Pfnt 



(18) 



R:,{xf) = {R'{xf,Pt)) . (19) 

In a first sub-section the production of anti-protons in these samples will be studied, as 
the determination of the net proton yields has to rely on the estimation of the cross sections of 
pair-produced protons. 

11.1.1 Anti-protons with final- state net baryon constraint 

The ratios RI^{xf) and RI^^{xf) for anti-protons are presented in Fig. [66] for forward 
(fT3l) and backward (fT5l) proton selection. 

For both net proton selections, a distinct correlation between the anti -proton yields and 
the presence of a tagged proton in the respective hemisphere is evident. The yield ratios reach 
values of 1.55 in the far backward and forward directions, respectively. The excess over 1 is 
halved dXxF = and ratio goes to 1 in the respective opposite xf regions. It has been verified 
that this correlation is not induced by asymmetries in the azimuthal acceptance of the detector. 
It is therefore due to an isospin effect as expected from the iso-triplet nature of baryon pair 
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Figure 66: Yield ratios Rf^i^xp) and RI^{xf) for anti-protons for a) forward and b) backward 
proton selection, as a function of xp- The full lines are forward-backward mirror symmetric and 
represent and estimate of the p overlap function 



production, Table [121 Indeed the presence of a net proton will unbalance the isospin structure 
present in the inclusive event sample towards the /s = -1 component of the corresponding heavy 
mesonic state and thereby enhance the p n yield in the selected hemisphere. From this argument 
follows a strict prediction for the tagging with net neutrons: in this case the I3 = +1 component, 
hence the np combination, should be favoured, and the an ti -proton yield should be reduced 
accordingly with respect to the inclusive sample. This is indeed the case as shown in Fig. [67] for 
forward neutron tagging, see ([141) . 




Figure 67: Yield ratios R^^^{xf) and RI^{xf) for anti-protons for forward neutron tagging, as 
a function of Xi;'. The full line represents the prediction using the extracted overlap function, 
Fig.|66l 
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This finding has two interesting consequences: it allows for the first time the estab- 
lishment of the range and shape of the anti-proton overlap function in a completely model- 
independent way. Secondly it has direct consequences for the possible production mechanism 
of baryon pairs in hadronic interactions. 

The overlap function is related to the full lines presented in Figs. [66]and[67l As the 
yield ratios are stable at \xf\ > 0.2 this means that the p feed-over is limited to \xf\ < 0.2. In 
the application to measured inclusive p yields the function should be normalized for projectile 
fragmentation to unity in the region xp > 0.2, should reach 0.5 at xp = and vanish at xp < - 
0.2, with anxp symmetric shape for target fragmentation. The resulting functions are presented 
inFig.lMl 
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Figure 68: Anti-proton overlap function as a function of for projectile fragmentation (full 
line) and target fragmentation (broken line) 

Concerning the production mechanism of baryon pairs, the strong isospin correlation 
both with the final state net baryons and with the hadronic projectile in the initial state speaks 
against the central production from the quark-gluon sea, as for instance gluon fusion. Concern- 
ing current hadronization models using string fragmentation, the ad-hoc introduction of diquark 
systems in the colliding baryons is necessary in order to describe the net baryon production. In 
addition quark/anti-quark and diquark/anti-diquark pickup processes have to be introduced for 
baryon pair production, in close resemblance to the description of hadronization in e+-i-e anni- 
hilation, with a multitude of adjustable parameters and doubtful predictive power, see Sect. ll 1.4l 

On the other hand the presence of heavy, high spin mesonic states in the early stage of 
hadronization of the highly excited baryonic systems in the p-i-p interaction presents a natural 
explanation for the observed correlations. In fact most observed heavy meson resonances have 
a p p decay branching fraction, or have been discovered in the inverse p-i-p annihilation process. 
It is interesting to note that baryon pair production, via the high effective mass involved, probes 
a rather primordial phase of hadronization as compared to the lighter final state hadrons. This 
again favours the observed strong isospin correlations. In this context it should be recalled that 
a high mass mesonic origin of baryon pairs was indeed proposed as early as 35 years ago by 
Bourquin and Gaillard |[26l in order to describe the observed inclusive yields. 

11.1.2 Pair-produced protons 

Following the isospin structure of baryon pair production, see Table [l2l anti-protons 
(p) and pair produced protons (p) are coupled in their yields by isospin symmetry with the /a 



pp^p 
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component of either the projectile or the trigger net baryon in the respective hemisphere. This 
is shown schematically in Fig. [69] where the p and p yields are presented with respect to the 
inclusive anti-proton level in p+p interactions. 
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Figure 69: Yield levels of anti -protons (p) and pair produced protons (p) with respect to the 
inclusive p yield in p+p interactions, left: for a change of projectile from +1/2 to -1/2, right: 
for a change of the trigger net baryon Is from +1/2 to -1/2 



Evidently the definition of "net" protons, which needs the knowledge of pair-produced 
proton yields, is linked to the measured p yields in a non-trivial fashion. As far as the overlap 
function between the forward and backward hemispheres is concerned, it is of course the same 
for p and p and has to be referred to the yield levels indicated in Fig. [691 



11.1.3 Net proton feed-over 

With the above clarification of the yields of pair produced baryons, the inclusive net 
proton yield may now be defined as: 

p'';\xf,pt) = Pr'(xi.,PT) - f^f{xF,PT) (20) 

As p-'^' is 1.6 times higher than p^^\ Fig. [69l this means a substantial decrease of the 
central net proton density with respect to the simple subtraction of the p yield. The consequence 
of this for the evolution of p^^^ with ^/s has been demonstrated in [12J where it has been shown 
that the central net invariant proton cross section approaches zero in the ISR energy range. 

The net proton feed-over and the corresponding overlap function may be determined, 
as for the anti-proton case Sect. II 1.1.11 above, by fixing a net proton in the projectile or target 
hemisphere following the selection criteria, ([T3]) and ([T5]) above, respectively. This ensures the 
absence of net protons in the corresponding hemisphere and results in the constrained net proton 
density: 

pI^'''\xf,pt) = p;{xf,Pt) - pI{xf,Pt) (21) 



These density distributions may be normalized by dividing by the inclusive yield, ([201) . 
resulting in the ratio 

DC,net/ „ ^ _ Pp'""(^F,Pr) 
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In a first instance, it can be shown that R^'^^^^xfjPt) is independent on pt in the range 
0.\<pt< 0.7 GeV/c, see Fig.|70] 
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Figure 70: Ratio R^'^^^{xf, pr) between constrained and inclusive net proton densities as a func- 
tion of Pt at different values oixp for net proton trigger in the projectile hemisphere (left side 
panels) and for trigger in the target hemisphere (right side panels) 



The xp dependence for the p^ integrated ratio i?p'"'^'(xir) and the px averaged ratio 
{R^^^^^xf.Pt)) are shown in Fig.|7T] 

Evidently the overlap functions for protons and anti-protons turn out to be the same 
within the statistical errors of the measurements. 

11.2 Proton and anti-proton feed-over in p+C interactions 

The study presented above for p-i-p collisions may be repeated for p-i-C interactions, 
albeit with reduced statistical significance due to the smaller available data sample. The density 
ratio for anti-protons, R^{xf) with forward proton tagging is presented in Fig. 1721 where the full 
line represents the feed-over function measured in p+p interactions, see Fig. [661 

Evidently the measurement reproduces the shape and extent of the p overlap in p+p 
collisions. 
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Figure 11: pr integrated and pt averaged constrained net proton density ratios _Rp'"''' as a func- 
tion of xp, a) for forward net proton selection and b) for backward net proton selection. The full 
lines shown represent the overlap functions presented in Fig.[ 




Figure 72: Anti-proton density ratio R^{xf) as a function of xp- The full line represents the 
feed-over in p-i-p events 



Net proton tagging in the backward hemisphere is not possible for p-i-A interactions 
as protons from intranuclear cascading ("grey" protons) which prevail in this area dilute the 
correlation. 

For the net proton density p^'^'^^xp) with tagging in the projectile hemisphere, a result 
which might look surprising at first view is found as shown in Fig. 1731 

Indeed the constrained net proton density reproduces the one found in p-i-p interactions 
down to xp = -0.2. This reveals an additional internal correlation effect generally present in 
minimum bias p-i-A interactions. By selecting a rather forward proton, see (fT3]) . single collisions 
corresponding to peripheral interactions are enhanced. In fact 60% of the minimum bias p-i-C 
events correspond to single interactions of the projectile proton with nuclear participants [|5]|. 
The situation is clarified by Fig. 1741 which presents the inclusive proton yield dn/dxp{xp) in 
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Figure 73: Net proton density in p+C collisions with forward net proton constraint as a function 
of Xp. The full line corresponds to the net density in p+p interactions 



p+C collisions (full line), 0.6 x dn/dxpixp) from p+p (broken line), and their difference (dotted 
line). It is evident that the bulk of protons in the tagging region (hatched area) comes from 
single collisions, with the fraction of multiple collisions in addition biased against very inelastic 
interactions which would favour \ow-xp protons by enhanced baryon number transfer. It may 
however be stated that also for p+C interactions with forward proton tagging the constrained 
net proton density dies out at Xi? > 0.2. 



Figure 74: The inclusive proton yield dn/ dxp^xp) in p+C collisions (full line), 0.6 x 
dn/dxpixp) from p+p (broken line), and their difference (dotted line). The inclusive proton 
yield dn/dxpixp) in p+p collisions is shown as well with full line. The hatched area represent 
the tagging region 



11.3 A remark concerning resonance decay 

The model-independent extraction of the baryonic feed-over and the experimental de- 
termination of the corresponding overlap functions, as described in the preceding sections, may 
be extended to resonance decay. As in fact most if not all final state baryons stem from res- 
onance decay it is of interest to investigate the consequences of the hadronic two-component 
mechanism observed for the final state baryons also for the parent generation of resonances. 
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This will be demonstrated below using as examples the well measured A++ resonance for the 
proton and a mesonic resonance approximated by measured tensor meson characteristics for the 
anti-proton feed-over. 



11.3.1 Protons from A.^^ decay 

The pt integrated density distribution dn/dxp^xp) of the A++(1232) resonance has 
been rather precisely measured by a number of experiments in the SPS and ISR energy ranges 
as shown by two examples [i27ll28ll in Fig. iTSk . The integration of the interpolated full line 
yields a total inclusive cross section of 7 mb. The corresponding decay proton distribution from 
— ^ p + 7r+ (100% branching fraction) is presented in Fig.lTSb together with the measured 
total inclusive proton yield as measured by NA49 Q, multiplied by the factor 0.27. 




Figure 75: pr integrated density distributions dn/dxp as a function of a) for A++(1232), b) 
for the decay protons from A++ — )• p + 7r+ (full lines). The broken line in panel b gives the total 
inclusive proton yield multiplied by a factor 0.27, the dotted line the difference between the two 
distributions. Also indicated in panel b is the 1/Mj distribution (see text) as dash-dotted line 



This plot demonstrates several important features: 

- the shape of the decay proton distribution from A++ follows, for < 0.6, very closely 
(to within about 5%) the total inclusive proton distribution re-normalized by a factor 
0.27. This means that 27% of all protons in this xp range stem from A++ decay alone. 

- the difference between the decay and inclusive distributions at > is well described 
by a form typical of single diffraction with ~ s(l — xp), dash-dotted line in 
panel b). 

- Taking account of the other A states A^(1232) and A°(1232) with branching fractions 
of 66% and 33% into protons, respectively, this means that more than 40% of all non- 
diffractive protons cascade down from A resonances 

- the shape difference between the dn/dxp distributions of the mother resonance and 
the daughter protons exemplifies the very effective baryon number transfer towards the 
central xp region in resonance decay even for low-Q resonances like the A. 

A number of further comments are due in this context. The absence of charge and flavour 
exchange in inelastic p-i-p interactions at SPS energy implies that the observed A resonances are 
not directly produced at least in the primordial phase of target and projectile excitation. They 
rather turn up as decay products of N* resonances. This fact has been experimentally proven in 
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a number of high precision studies of single and double diffraction into ptt+tt states from PS 
to ISR energies Il29l - [32ll . In these final states which cover an la;^! range from 1 to about 0.6 Il32]| . 
the initial /, state is 1/2, +1/2. Although the sub-channels p7r+ and pvr^ show clear A(1232) 
signals, those are completely contained in the decay mass spectra of a series of N* resonances 
(N*(1440), N*(1520), N*(1680)). In the more central area of hadronization, this clear isospin 
signature will become diluted and contain 1/2, -1/2 states, see also the isospin correlations 
discussed above. Nevertheless the production and decay of N* resonances without intermediate 
A states will provide another important source of final state nucleons. 

The two-component mechanism of hadronization should therefore also be considered 
for baryon resonances with the constraint to reproduce the measured overlap function for final 
state protons. This is fulfilled for A^"*" by the separation of the production cross section into a 
target and a projectile component as indicated in Fig. r76k. The corresponding xp distributions 
for the decay protons are given in Fig.r76b. 




Figure 76: Density distributions dn/dxp as a function of a) for the target and projectile 
components of A++ and b) for the decay protons 



This choice results in the ratios R = {dn/dxF)pToi/ {dn/dxp)"^^^ between projectile com- 
ponent and inclusive distribution shown in Fig. 1771 




Figure 77: Density ratios R = {dn/ dxF)pm]/ {dn/dxpY^^^ a) for A"*""*", full circles and hatched 
line, and b) for the decay protons from A"*""*" decay, open circles, as a function of xp- The 
measured overlap function, see Fig.jTTl is indicated as the full line in both panels 
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Evidently the measured proton feed-over is precisely reproduced, Fig.lTTb. if the overlap 
function of the mother resonance is chosen slightly wider in xp. Fig. ITTk . indicating a certain 
mass dependence in the xp scale. In this context the dependence on decay particle mass exhib- 
ited by the much reduced feed-over for the decay pions from as elaborated in [5}, should 
be recalled here. 



11.3.2 Anti-protons from heavy meson decay 

A complementary approach may be used for the production of anti-protons from heavy 
mesonic states |l26l|. In fact a sizeable number of states above the pp mass threshold at 
1.88 GeV/c^ have been observed ll33l both in hadronic interactions and in the inverse pp anni- 
hilation into final state hadrons. If due to the large width and density of these states their direct 
experimental detection in mass spectra is difficult, the application of Partial Wave Analysis per- 
mits their localization and determination of quantum numbers. In the following a hypothetical 
state at 2.5 GeV mass with a Breit-Wigner width of 0.25 GeV and a two-body decay into pp 
is used in order to study the corresponding daughter xp distributions and overlap functions. 
An invariant xp distribution consistent with the one for vector and tensor mesons published by 
Suzuki et al. [34] yields the inclusive density distribution dn/dxp for the decay anti-protons 
shown in Fig. 1781 




Figure 78: a) Invariant xp distribution F{xp) of a mesonic state with 2.5 GeV mass (full line) 
compared to various vector and tensor mesons [34 1 normalized to unity dX xp = 0,h) resulting 
inclusive density distribution dn/dxp for the decay anti-protons (full line) compared to the 
measured yields in p+p interactions [|3]| (broken line) 



Evidently this choice reproduces perfectly the measured inclusive anti-proton yields 
measured by NA49 [3]. Imposing the target-projectile decomposition for the yield distribution 
of the heavy meson as shown in Fig.lTQb. the overlap function for the decay baryons reproduces 
closely the measured feed-over for anti-protons, Fig.lTOb. 

As far as the normalization of the resulting inclusive yields is concerned, it may be stated 
that - compared to the p cross section of 1 .2 mb |[3l - the measurement of [34] indicates 1 .6 mb 
for the f4(2050) state alone, and the mass dependence of the tensor meson production given 
in fSS] an f4(2050) cross section in excess of 5 mb at ISR energy. There are, however, almost no 
measurements of the branching fraction into baryon pairs, and if so, they vary by large factors, 
for the f4(2050) for instance from 50% [[36l to 0.2% [[37ll . The amount to which anti-protons 
cascade from heavy meson decay has therefore to stay an open question. 
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Figure 79: a) Density distribution dn/dxp for the heavy mesonic state (full line) normalized 
to 1 at a; i? = decomposed into a target and a projectile component (broken lines), b) the 
corresponding overlap function R(xf) (full line) superimposed with the measurement (open 
circles and broken line) 



11.4 Comparison to a microscopic simulation code 

In the preceding NA49 publications concerning p+p and p+C interactions [[BQ a strictly 
model-independent approach to the interpretation of these extensive and precise data sets has 
been followed. Comparison to the multitude of existing microscopic simulation codes has there- 
fore been avoided in clear appreciation of the fact that the actual understanding (or, better, lack 
of understanding) of the theoretical foundations of the soft sector of QCD calls for improved 
experimental information rather than ad hoc parametrizations. Given the detailed data concern- 
ing baryons both on the inclusive level and on the level of baryonic correlations contained in ||3l 
and in the above discussion, it might however be useful to confront the obtained results with 
one specific microscopic simulation in order to obtain some idea of the precision and predictive 
power reached in such approaches. 

The relatively recent code named JAM [l38l| has been chosen for this comparison. This 
code is supposed to describe hadronization in elementary and nuclear interactions over the com- 
plete cms energy range from 1 GeV up to collider energies. It uses, in the SPS energy range, soft 
string excitation following the HIJING approach which in turn relies on the string fragmentation 
mechanism developed in the PYTHIA/LUND environment. 

11.4.1 Inclusive baryon and anti-baryon density distributions dn/dxp 

The measured proton and neutron density distributions dn/dxp in p-i-p interactions are 
confronted in Fig. [80] with the predictions from JAM. 

Large systematic deviations between prediction and data are visible at low xp (-40%), 
at medium xp ^ 0.25 (-1-50%) and at high xp (-30%) for the protons. Systematic deviations of 
similar or bigger size are also seen for the HSD and UrQMD codes [|39l . 

The neutron yields are evidently obtained from the proton densities by a constant multi- 
plicative factor of 0.69 with the exception of the large xp region where a diffractive component 
with a behaviour, see Sect. II 1.3.1[ is added to the protons. The deviations are -1-60% at 

medium xp and more than -100% at large xp for the neutrons. A look at the density distri- 
butions for n-i-p interactions. Fig. [STl shows only approximate isospin symmetry which would 
impose that neutrons from proton fragmentation should be equal to protons from neutrons in the 
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Figure 80: Density distributions dn/dxp as functions of a;i? a) for protons and b) for neutrons 
(full circles) in p+p collisions compared to the prediction from JAM (full lines) 



xp regions beyond the target-projectile overlap. The same should of course be true for neutrons 
from neutron beam and protons from proton beam. 




0.5 1 0.5 



Figure 81: Density distributions dn/dxp as functions of xp a) for neutrons and b) for protons 
in n+p collisions from JAM 



The anti-baryon densities shown in Fig. [82] show an interesting pattern, the anti-proton 
and anti-neutron yields being identical for both projectile-target combinations. 




Figure 82: a) anti-proton and b) anti-neutron distributions dn/dxp as functions of xp from p+p 
and n-i-p interactions. The anti-proton distribution measured by NA49 is shown in panel a) with 
full circles 
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This means that only the /s = combinations pp and nn are allowed which are produced 
with equal yields. The comparison with the measured anti-proton density in p+p collisions 
shows a sizeable underestimation of the yield by 35% at xp = and an equally large overesti- 
mation atxp ~ 0.3. Invoking the isospin effect measured in n+p interactions [i l2il this difference 
will increase to 50% in this reaction. 

11.4.2 Baryonic correlations 

If it might be a straightforward possibility to remedy the observed discrepancies between 
the prediction and the inclusive data by modifying some of the many parameters involved in the 
simulation codes, the correlation data will probe the "physics" input on a deeper level. This 
applies especially to the measured isospin effects. 

A first comparison concerns the net proton density correlated with a trigger baryon in 
the projectile hemisphere resulting in the overlap function R^'^^^, (f22l) . as shown in Fig. [83] for 
the two trigger xp bins defined in Sect. ll l.ll Eqs.[T3]and[T4l 




Figure 83: Net proton overlap function R^'^'^^xp) as a function of xp for a) proton trigger at 
xp = 0.35 to 0.5, b) neutron trigger atxp = 0.35 to 0.5, c) proton trigger atxp = 0.5 to 0.7 and 
d) neutron trigger at Xi;' = 0.5 to 0.7. The full lines represent the measured function. Sect. II 1.2[ 
the points come from the JAM simulation code 

The microscopic simulation results in a feed-over behaviour which only reaches to about 
±0.06 in Xp. This is considerably sharper than the measured behaviour (full lines in Fig. [83]) 
and corresponds to the pionic feed-over extracted in 10. In addition there is an asymmetric long 
range tail that extends up to and beyond I^fI ~ 0.3 and is different both for proton and neutron 
trigger and for the two trigger x p bins in contrast to the data. 

The anti-proton feed-over behaviour characterized by the ratio R^{xp) of correlated to 
inclusive densities, Eq[T8l is presented in Fig.[84]again for the two available trigger Xi? bins and 
for proton and neutron trigger. 

A rather complicated pattern concerning the simulation emerges. Evidently the sym- 
metric isospin effect observed experimentally between proton and neutron triggers is not re- 
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Figure 84: Anti-proton density ratio Ry^xp) as a function oixp for a) proton trigger dXxp = 0.35 
to 0.5, b) neutron trigger at xp = 0.35 to 0.5, c) proton trigger at xp = 0.5 to 0.7 and d) neutron 
trigger at = 0.5 to 0.7. The full lines represent the measured function, Sect. lll.lTTl the points 
come from the JAM simulation code 



produced although there is a general reduction of the density ratio in the trigger hemisphere 
with neutron triggers. As only isospin singlet baryon pairs are generated, see Sect. II 1.4.1[ any 
isospin effect is not really expected. The strong suppression of the density ratio in the high-Xi? 
trigger bin already starting at xj? ~ 0.2 indicates an effect of energy-momentum conservation 
in the baryon pair simulation which is probably the result of the details of string fragmentation. 

In conclusion the microscopic simulation results in major deviations from the data both 
on the level of the inclusive and of the correlated yields. The absence of a proper treatment of 
isospin effects both concerning the component of the projectile and of the final state baryons 
is flagrant. This puts into doubt the application of this approach to p-i-A and especially A-i-A 
interactions where the neutron component is preponderant. 

12 Hadronization in p+C collisions: anti-protons 

The two-component hadronization mechanism studied experimentally in the preceding 
section allows, in conjunction with the precision data on anti-proton production in p-i-p inter- 
actions published by NA49 [3|, the confrontation of the measured p yields in p-i-C collisions 
with a straight-forward prediction based on elementary reactions. This is simplified in the case 
of anti-protons by the fact that there is, in contrast to proton and pion production, no contribu- 
tion from nuclear cascading (see Sect. |9] above). In consequence the superposition of target and 
projectile fragmentation should suffice to completely describe the observed cross sections. In a 
first step it may be assumed that the projectile contribution corresponds exactly to the one in 
p+p interactions. For the target contribution the same basic assumption may be made with two 
important additional constraints taking account, firstly, of the multiple intranuclear collisions of 
the projectile, and secondly, of the isospin factor involved in the fragmentation of the neutrons 
contained in the nucleus. 
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The mean number of projectile collisions {v) in minimum bias p+C collisions has been 
investigated in [|3 using three different methods: 

- a Monte Carlo simulation using the measured nuclear density profile as input with the 
result (z/) = 1.6 

- the measured p+C inelastic cross section which gives, via the relation 



A ■ (y{pp) 



(23) 



a{pA) 

an estimation of (z/) under the assumption that the intranuclear inelastic interaction cross 
section stays constant for all subsequent projectile collisions, resulting in (u) = 1.68 
- the measured increase of pion yields in the backward hemisphere using the fact that the 
contribution from intranuclear cascading as well as the one from projectile fragmenta- 
tion die out atxp ~ -0.1, yielding (z/) = 1.6. 

The isospin factor for p production from the isoscalar C nucleus may be calculated, 
using the measured increase of p yields in n-i-p collisions [lT2ll as 



0.5(1 + (7^) = 1.3, 



(24) 



with = 1.6. The combined overall factor to be applied to the target component is 1.6 x 1.3 = 
2.08. 



12. 1 pt integrated density dn / dxp (xp ) 

The evolution of the pt integrated anti-proton density dn/dxpixp) from the elemen- 
tary p-i-p to the p+C interactions using the superposition of target and projectile components 
discussed above is presented in Fig. 




Figure 85: Measured integrated anti-proton density dn/dxpixp) in p+C interactions (open 
circles) confronted with the superposition of target and projectile components in p+p collisions 
(dotted lines) and total yield (dash-dotted line), with multiplication factors {u) = 1.6 and 1 for 
target and projectile, respectively, broken line, and with the additional isospin factor 1.3 for the 
target component (full line) 



This most straight-forward superposition picture evidently reproduces the measured p 
densities quite closely, overestimating them by about 10% in the target area and underestimating 
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them by the same amount in the projectile region. This can be remedied by decreasing the target 
multiplication factor by 12% from 2.08 to 1.84 and by increasing the projectile contribution by 
10% resulting in a reasonable fit of the experimental data compatible with their statistical errors 
as shown in Fig. [861 




Figure 86: Measured pr integrated anti-proton density dn/ dxp^xp) in p+C interactions com- 
pared with a target-projectile superposition scheme with factors 1.84 and 1.1, respectively, for 
the target and projectile components. The yield from p-i-p is shown with dash-dotted line 

The increase of the p yield by 10% in the projectile hemisphere with respect to the direct 
estimation from p-i-p collisions is a first interesting consequence of this study. It is equal to the 
increase of the pion yields deduced in [5J. 

The reduction of the target contribution by about 12% with respect to the simple super- 
position of elementary hadronization processes as characterized by the mean number of colli- 
sions (z/) and isospin symmetry, is a second important result. The magnitude of this reduction 
being relatively small it is nevertheless on the limit allowed by the experimental determination 
both of (z/) [j5l and of the isospin effect on anti-proton production from neutrons [12]. 

Regarding the excitation mechanism of colliding hadrons by the exchange of gluons or 
gluonic (charge and flavour-less) objects, the observed effects are however to be expected. The 
projectile interacts in its multiple collisions subsequently with "fresh" nucleons which did not 
undergo previous exchanges. Hence its effective excitation level will increase with {u). The 
projectile on the other hand suffers in each collision a loss of its gluonic component such that 
less excitation energy with the target nucleons can be exchanged in subsequent interactions. In 
this sense the nucleus may be regarded as a gluon filter, the description of this phenomenology 
by the term of "energy loss" giving only a very general and somewhat misleading impression. 

The anti-proton yields regarded here are especially sensitive to this notion as at SPS 
energy the s-dependence of the production cross section is still rather steep with about 10% per 
GeV in ^/s [|26l . This is in contrast to the production of mesons, with 5% per GeV for mean 
kaons ^ and only 3% per GeV for mean pions [1401 . In this sense the study of pion yields in 
the target hemisphere approaches a precise measure of {u) [[51 whereas the observed anti-proton 
cross section indicates an effective loss in -y/s of about 1 GeV in target fragmentation for p-i-C 
collisions. Evidently the extension of this study to heavier nuclei is of considerable interest in 
this respect. The data on p-i-Pb collisions with controlled centrality available from NA49 will 
illuminate this point, as is already visible in the preliminary results shown in |[4T1l . 
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12.2 Double differential invariant densities f (xp , Pt) / o^inei 



Further details of the two-component hadronization mechanism detailed above become 
visible if comparing the double differential invariant densities to the optimized prediction from 
target and projectile fragmentation obtained from the px integrated yields, see Fig. [861 This is 
shown in Fig. [87] where the invariant cross sections f{xp,pT) per inelastic event (open circles) 
are presented together with the predicted densities (full line). 




-0.2 0.2 




Figure 87: Double differential invariant anti-proton densities per inelastic event (open circles) 
compared to the prediction from the two component fragmentation mechanism (full lines) as 
a function of xp, for different values of pr between 0.1 and 1.3 GeV/c. The distributions for 
different pt values are successively scaled down by 2 for better separation 
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Three main features may be extracted from these plots: 

- the target component reproduces, within the experimental errors, the densities predicted 
from elementary interactions forpT < 1 GeV/c. 

- the projectile component reveals a definite px dependence. The measurements fall below 
the prediction for p^- < 0.5 GeV/c and increase smoothly above the prediction above this 
value. 

- in the range of > 1 GeV/c the excess of the projectile component starts to extend 
well into the target hemisphere down to Xi;' ~ -0.2. 

This behaviour is detailed in Fig. [88] which presents the ratio R between measurement 
and prediction for the px values shown in Fig. [87j 
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Figure 88: Ratio R between measured and predicted anti-proton densities for different pt values 
between 0.1 and 1.3 GeV/c as a function of xp- The full lines represent the parametrization 

R={l-F^) + RP'°iF^ 



The same ratio is plotted against pt for different values of xj? in Fig. [89l 
It should be stressed here that these experimental results are in strong support of the 
independent target-projectile fragmentation in p-i-C interactions. The modification of the pr 
distribution of the projectile component which superimposes itself to the 10% increase in total 
yield, indicates that the Cronin effect, whose onset is visible here, is limited to the projectile 
hadronization. The increase of yields at higher p-p which extends well into the target hemisphere 
might be due to an extension of the feed-over range with pt (see also (H for pions), although the 
limited statistics does not allow for quantitative statements. The increase of R in the projectile 
hemisphere may be parametrized as 

i?p>-°j = 0.9 + 0.35p| (25) 

Its modification in the transition to the target hemisphere is then predicted by the pro- 
jectile overlap function (Fig. I68l) resulting in the full lines indicated in Figs. [88] and [891 

Further information is contained in the detailed comparison of the average pr values to 
p-i-p interactions given in Fig. [90l 
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Figure 89: Ratio R between measured and predicted anti -proton densities for different x p values 
as a function of px- The full lines represent the parametrization R = {1 — F-^) + R^^°^F^ 
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Figure 90: a) comparison of average px between p-i-p collisions (full line) and p+C interactions 
(data points), b) difference in (px) between p-i-C and p-i-p interactions. The full line in panel b) 
represent the anti-proton overlap function F-^ from Sect. II 1.1.11 



The convergence of the difference in (px) towards zero in the backward hemisphere, 
Fig.l90b. and its description by the overlap function F-^, Fig. [681 is to be regarded as yet another 
manifestation of the two-component mechanism of hadronization as discussed in Sect. [TT] 

13 Proton production in p+C collisions: pr integrated yields 

As in the preceding section on anti-protons, proton production will be first discussed 
using the px integrated yields dn/dxp in order to clearly visualize and separate the three basic 
contributions to the overall proton cross section. After establishing the net proton density by 
subtracting the yield of pair produced protons, the fragmentation of the hit target nucleons, the 
projectile fragmentation and the contribution from intra-nuclear cascading will be treated in 
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turn. 



13.1 Pair produced protons and net proton density 

It has been shown above that the iso-triplet nature of baryon pair production imposes a 
detailed follow-up of isospin effects all through the hadronization process. This means of course 
that the yield of anti-protons is in general not identical to the yield of pair-produced protons. It is 
therefore not sufficient to subtract the anti-proton density as established in the preceding section 
from the proton yield in order to obtain the net proton density. In the case of the asymmetric 
p+C interaction, the equality of anti-proton and pair-produced proton yields is only valid for 
the target part of the overall cross section due to isoscalar carbon nucleus, invoking isospin 
symmetry. For the projectile component on the other hand this does not apply. In the absence 
of precision data for anti-proton production in n-i-C collisions it may be assumed that the same 
enhancement of pair-produced protons over anti-protons as in p-i-p interactions applies where a 
factor of 1.6 has been established lfT2l . This modifies the superposition scheme for anti-protons 
described in Sect. [I2lby increasing the projectile contribution from 1.1 to 1.76 with respect 
to the input p+p densities. The resulting px integrated density distribution for pair produced 
protons (hereafter denoted as p ) is shown in Fig.l9Tl 
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Figure 91: Density distribution dn/dxpixp) for pair produced protons (p ) as a function of 
(full line). The data points are the measured anti-proton densities 

This distribution is in contrast to Fig. [80l almost symmetric and the comparison to the 
measured anti-proton densities demonstrates that this isospin effect is by no means negligible. 

The corresponding net proton density distribution results from the subtraction of the pair 
produced protons from the total proton density distribution. Table |7] and Fig. |52l It is presented 
in Fig. |92]in the x p range influenced by pair production, -0.4 < xp < +0.4. 

The net proton density shows, in contrast to the total proton yield, a smooth behaviour 
around xp = 0. 

13.2 Target and projectile components in net proton production 

The fragmentation of those nucleons in the carbon nucleus which are hit by the through- 
going projectile, here called the "target component", is a quantity which should be closely 
related to the proper superposition of net proton production in the elementary p-i-p and p-i-n 
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-0.4 -0.2 0.2 0.4 

Figure 92: Density distribution dn/dxF{xp) for protons (broken line) and for net protons (full 
line) as a function of xj? in the range -0.4 < xp < +0.4 

interactions, multiplied by the number of projectile interactions. Such a superposition has been 
shown to describe the target component of the anti-proton yields, Sect. [121 up to a loss of 
about 12% specific to multiple collisions and related to the strong s-dependence of baryon pair 
production. For baryons this loss should be negligible as the baryon density is to first order 
s -independent at SPS energy. 

For the prediction of the target component, knowledge about the two basic contributions 

from 

p-i-p — p and p-i-n — )■ p (26) 
is needed. The latter process may be related, via isospin symmetry, to the reaction 

p-i-p n (27) 

which is measured by the NA49 experiment [3|. The resulting prediction for the isospin aver- 
aged net proton density from p-i-p collisions is described in the following section. 

13.2.1 Isospin averaged net proton density from p+p collisions 

The pt integrated proton and neutron densities dn/dxp as measured by the NA49 ex- 
periment are shown in Fig. |93]as dotted and broken lines. The corresponding net proton yield is 
obtained by subtracting 1 .6 times the measured anti-proton yield. For the net neutron yield, the 
subtraction of the measured anti-proton yield is indicated by the isospin symmetry of baryon 
pair production. In the case of the isoscalar carbon nucleus, a simple average of the two distri- 
butions is to be performed as indicated by the full line in Fig.l93l 

13.2.2 Predicted target component of net proton production 

In order to obtain the predicted target component of net proton production in p-i-C inter- 
action, the predicted yield from p-i-p collisions has to be multiplied by the number of projectile 
collisions (u) = 1.6 and by the target overlap function. Sect. II 1.21 The resulting distribution is 
shown in Fig. [941 together with the total net proton yield. Sect. 113. 11 

As a first result it may be seen that the target component nearly saturates the total yield 
at xp ~ -0.2. 
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Figure 93 : Total and net proton density dn/ dxp^xp), dotted lines , total and net neutron densities 
(broken lines) together with the measured neutron cross sections (data points from |I3|) and the 
isospin averaged net proton density(full line), as a function of 




Figure 94: Total net proton density dn/dxpixp), broken line, and predicted target contribution 
(full line) as a function of 



13.3 The projectile component of net proton production 

The subtraction of the predicted target component from the overall net proton density. 
Fig. mi allows now for the extraction of the projectile component of the p+C interaction. This 
is demonstrated in Fig. |95] which shows the net projectile component (full line) in comparison 
with the net projectile component in p+p collisions (broken line. Sect. 111.2b . 

Several features of these distributions are noteworthy: 

- this study allows for the first time the isolation of the net proton projectile component in 
p+A collisions over the full phase space 

- in principle this net proton component is strongly constrained by baryon number con- 
servation in comparison to p+p collisions, to the extent that the surface under the two 
distributions should be equal up to second order effects like a relative increase of neutron 
or hyperon production. 

- in fact the integrated yields are 0.547 and 0.52 net protons per inelastic event for p+p 
and p+C interactions, respectively. In view of the multi-step methodology involved in 
extracting this experimental result, the agreement to within about 5% is certainly com- 



74 



"O 
"D 

0.6 
0.4 
0.2 


0.5 1 

Figure 95: Net proton density dn/dxpixp) for the projectile component in p+C interactions 
(full line) and in p+p collisions (broken line) as a function of xp- The shaded area indicates the 
onset of the contribution from the nuclear component in p+C interactions 

patible with the systematic uncertainties. 

- a relative increase of the net neutron yield, which is in principle not excluded in multiple 
hadronic interactions, is improbable as it has been shown, see Sect. ll0.3l that thebaryon 
number transfer is identical for neutrons and protons. 

- a relative increase of hyperon production can as well only have limited influence as an 
increase of A and S production by 50% would only reduce the observed difference from 
5% to 4%. 

- the shapes of the two distributions are rather similar, with a downward shift of about 
0.15 units in xp in p+C except for the diffractive region which is governed by single 
projectile collisions 

- at xp < -0.2 there is a steep increase of the target- subtracted density (shaded region in 
Fig.|95]) which is due to the tail of protons from intra-nuclear cascading. This contribu- 
tion will be discussed in detail below (Sect. 113.41) . It is however to be noted that a clean 
separation of the projectile component from both the target and the nuclear cascading 
contributions is being achieved at SPS energy. At lower cms energies, the nuclear com- 
ponent will extend into higher Xi;' ranges, covering xp = at AGS energy, as it scales in 
Piab rather than s^;^. 

It is interesting to extract the projectile overlap function from the ratio R^'^^\ see (l22l) . 
between the projectile component and the symmetrized total density. This is shown in Fig. |96] 
in comparison with the corresponding function in p+p collisions. 

As a further important result of this study it may be stated that in multiple collisions the 
width of the proton feed-over from the projectile to the target hemisphere widens. In view of the 
discussion of resonance decay in Sect. II 1.31 above, this would be compatible with an increase 
of the effective mass in the process of projectile excitation. 

13.4 The nuclear component of net proton production 

If the subtraction of the target component from the total net proton density distribution 
reveals the projectile component at a;^ > -0.2, it also should allow for the extraction of the 
nuclear component at xp < -0.2. Here the nuclear component is defined as the retarded proton 
density produced by the intra-nuclear cascading of nucleons following the momentum transfer 
from the projectile to the nucleus, as opposed to the prompt fragmentation of those nucleons 
which are hit by the projectile. 
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Figure 96: Projectile overlap function R^'^'^^xp) as a function of xp from p+C interactions (full 
line) compared to the same function from p+p interactions (broken line) 



The situation is clarified in Fig. l97l which shows the total net proton density dn/dxp 
for Xp < -0.2 (full line) together with the target component discussed in the preceding section 
(broken line), Fig.|94l 




Figure 97: Net proton density dn/dxp as a function of xp in the range -1.6 < < -0.2. Full 
line: total measured yield, broken line: predicted target component 



The subtraction of the target component from the total yield. Fig. |98] (full line), results 
in a rather particular pattern. 

This pattern features two distinct regions in xp. For the region of ±0.2 units around 
Xp = - I, an abrupt bipolar instability arises from the presence of a diffractive peak in the pre- 
dicted target component, whereas for the regions -1.6 < xp < -1.2 and -0.8 < xp < -0.2 a 
smooth Xp distribution results which is well fitted by a Gaussian of the form 

dn (xp+0.92)^ 

-. 0.95e 2.0.265^ (28) 



dx 
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Figure 98: Net proton density dn/dxp as a function of xp in the region -1.6 < xp < -0.2 
resulting from the subtraction of the predicted target density from the total measured proton 
yield (full line) and a Gaussian fit to the difference in the two regions -1.6 < xp < -1-2 and 
-0.8 < < -0.2, broken line 

centered at = -0.92 with an rms of 0.265 as shown by the broken line in Fig.[98l 

The diffractive component of the predicted target density should show up in any mea- 
surement of the total proton density in the region around Xi? = -1 if the nuclear component 
would have as expected a smooth behaviour across this region. Assuming the Gaussian fit to 
describe this smooth behaviour one may tentatively subtract the Gaussian shape from the total 
proton density as shown in Fig.l99l 




Figure 99: Net proton densities dn/dxp in the region -1.6 < xp < -0.2 as a function of xp, 
full line total measured yield, broken line Gaussian fit to the target density in the regions - 
1.6<xp'<-1.2 and -0.8 < xp < -0.2, dotted line difference of the two distributions 
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This results in a modified target component (dotted line in Fig. |99]) which reaches down 
to xf values at ~ -1.2 and goes smoothly through xp = -l. The difference between the modified 
and predicted target component is bipolar around xp = - I and conserves, to first order, the total 
integrated density, as shown in Fig. llOOl 
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Figure 100: a) Difference between the modified target component (Fig. |971 broken line) and 
the predicted component obtained by subtraction of the Gaussian nuclear component (Fig. [991 
dotted line), b) Gaussian longitudinal proton density distribution as obtained in quasi-elastic 
scattering of alpha-alpha collisions at a/s = 30 GeV/nucleon at the CERN ISR [|45ll 



In this context the action of the Fermi motion of the hit nucleons on the momentum 
distribution of protons in the diffractive peak should be recalled. Fermi motion may be approx- 
imated by a Gaussian momentum distribution of about 75 MeV/c rms width in the nuclear rest 
system. This translates at SPS energy to an rms width of about 0.085 in xp- Low momentum 
transfer processes like quasi-elastic scattering or diffraction convolute with this momentum dis- 
tribution. This leads to a smearing of the peak structure characterized by a base width of ±2-3 
times the Fermi rms around xp = -l. This is quantified by the Gaussian x p distribution with rms 
0.085 in Fig. 1 100b which covers exactly the range of the observed deviation from the predicted 
diffractive peak. This convolution has been studied in quasi-elastic alpha-alpha scattering at the 
CERN ISR 03| where it was shown that low momentum transfer re-scattering processes have 
only small influence on the longitudinal width in contrast to a substantial broadening of the 
transverse momentum distribution, see Sect. 14 below. 
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The extracted width of the nuclear component in longitudinal momentum with an rms 
of 0.265 units in a; is much larger than the spread expected from Fermi motion. This feature 
is known from quite a number of experiments since some decades [6] but is shown here for the 
first time to be symmetric, however not around xp = - I but around xp = -0.92. This substantial 
transfer of baryon number needs comment. One effective mechanism contributing to the transfer 
of nucleons in longitudinal momentum is certainly the production of pions in the intra-nuclear 
cascade. In the excitation of / = 3/2 states at pion threshold, the decay products are located at 
Xp = -0.85 and x^? = -0.15, respectively, for protons and pions, and for isobars at xp ~ -1. For 
isobar masses above threshold the xp values of the decay baryons will shift further up. 

The production of pions in the nuclear cascade has indeed been studied in the accompa- 
nying publication The backward pion yields have been separated into the two components 
of target fragmentation which may be predicted from elementary hadronic collisions, and of 
nuclear cascading, in close similarity to the argumentation carried out here for baryons. As a re- 
sult the total number of pions in the nuclear component has been determined to 0.3 per inelastic 
event. 

The Gaussian distribution of the nuclear component. Fig. |99l integrates to 0.6 protons 
per inelastic event, see Sect. ll3.6l below. Isospin symmetry and the absence of charge exchange 
in the primary p-i-C interaction at SPS energy ask for an equal number of neutrons to be pro- 
duced. This means that every fourth proton in the nuclear cascade produces a pion. In the decay 
of the isospin 3/2 isobars involved with pion production, the decay protons are transferred to 
Xp = -0.85 at pion threshold with the decay pions being centred at a; = -0.15 for low lab mo- 
menta. Integrating over the pion transverse momentum distribution this center shifts down to 
about -0.2 [16| with a rather long tail towards lower xp. The corresponding nucleons will be 
placed at Xp = -0.8 and above. The center of the nuclear component at xp = -0.92 is there- 
fore compatible with pion production in the nuclear cascade via isobar excitation near threshold 
and the subsequent decay into baryons and pions. The symmetry of the nuclear component 
would however ask for a bigger pion yield at about one pion for every second proton. Here 
de-excitation processes without pion emission like final state isobar-nucleon re-scattering could 
be involved to reduce the pion yield. 

In conclusion of the preceding sections on pt integrated proton distributions it may be 
stated that a three-component mechanism of convincing internal consistency has been estab- 
lished in a quantitative and precise way. These three components are defined by 

- the fragmentation of the projectile which obeys baryon number conservation and shows 
a transfer in longitudinal momentum corresponding to 0.15 units in Feynman xp. 

- the prompt fragmentation of the nucleons hit by the projectile. This component has been 
predicted from elementary proton-nucleon interactions invoking isospin symmetry and 
the mean number of projectile collisions. 

- the nuclear component which arises from intra-nuclear cascading and is partially ac- 
companied by pion production. 

At SPS energy the nuclear component extends up to xp = -0.2. It is well separated 
from the projectile fragmentation which feeds over into the target hemisphere down to the same 
value of Xp = -0.2. The target fragmentation in turn reaches from xp = +0.2 down to xp = - 
1.2. The latter range is well beyond the kinematic limit of elementary fragmentation with a 
diffractive peak close to and above xp = - I. This sharp diffractive structure is evidently diluted 
and smoothed by quasi-elastic re- scattering of its low-momentum nucleons in accordance with 
experimental results from other experiments. 
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13.5 A remark about s-dependence 

If the SPS energy range is, as shown above, a well suited region for the separation of the 
three components contributing to the measured net baryon density, it is interesting to look at the 
expected evolution with y/s. 

13.5.1 The nuclear component 

The study [6] of backward proton production in p+C interactions shows that only small 
changes in the momentum distributions are to be expected by increasing -y/s to RHIC energy 
and beyond, or by decreasing it into the AGS energy region. Indeed one has to move down to 
^/s below about 3 GeV in order to see threshold effects drastically reducing the proton density. 
In terms of the range mxp covered by the nuclear component, however, kinematics will extend 
the upper limit of ~ -0.25 at SPS energy to higher values until a;^? = is reached at ~ 4 GeV. 
This purely kinematic effect is presented in Fig. llOll where a^i? is plotted as a function of l/-\/i 
for two values of piab and three values of pr- 



X 



0.5 


Pt 


1 1 1 1 1 
= 0.0 GeV/c 






Pt 
-Pt 


= 0.2 GeV/c 
= 0.4 GeV/c 


p =2GeV/c / 

lab // 























,_,;;;======>^^p = 1 GeV/c ' 










0.5 






1 < , < < 1 < < < 



0.1 0.2 0.3 0.4 

1/Ms[GeV-^] 

Figure 101 : a; as a function of l/-\/i for p\^^ = 1 and 2 GeV/c, varying from to 0.4 GeV/c 

This behaviour means that for given xp (or rapidity) the nuclear contribution will in- 
crease with decreasing beam momentum. This effect is of course also to be expected for pe- 
ripheral heavy ion collisions where the separation of prompt baryons from the delayed nuclear 
cascade will become more important but also more difficult with decreasing interaction energy. 

13.5.2 The target fragmentation 

As shown above prompt net baryons from target fragmentation are well described by 
the superposition of elementary proton-nucleon collisions. With increasing -y/i this means that 
this contribution will move back in xf due to increasing transparency until at the highest ISR 
and at RHIC energies the central region around xp = will be only populated by pair pro- 
duced baryons. In the xp region below ~ -0.3 the yields will not change appreciably due to 
the approximate scaling behaviour. Towards lower -\/s the feed-over of the target fragmentation 
into the projectile hemisphere will accordingly increase. Together with the increasing nuclear 
component a clear separation will therefore become more involved. 
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13.5.3 The projectile component 

The evolution of the projectile component with collision energy is of considerable inter- 
est as it carries unique information about multiple hadronic collisions. Since its diffractive part 
is connected - for minimum bias p+A collisions - to the fraction of single projectile collisions 
which is as discussed above on a level of 60% for p+C interactions, there should be approxima- 
tive energy scaling with a progressive widening of the peak towards lower energies as measured 
in p-i-p interactions. For the non-diffractive part the situation is however much less clear, in par- 
ticular as there is no theoretical understanding of the observed baryon number transfer and as 
there is only very limited experimental information available on the required level of precision. 
It has been shown above that the net baryon density distribution shifts downward by 0.15 units 
of in p-i-C collisions, the same shift being observed for the neutrons. It has also been shown 
that this baryon number transfer increases with centrality as indicated by the dependence on the 
number of grey protons. It has also in addition been indicated that resonance production and 
decay presents a very effective source of baryon transfer. If with increasing interaction energy 
the excitation of the projectile proton in its multiple collisions also increases and therefore the 
spectrum of produced resonances extends to higher masses, the reduction of central net baryon 
density via transparency should be reduced or even compensated. At lower y/s the situation 
should become, as for the two other contributions, more obscure due to the increased mutual 
overlap of the different components. 

13.6 Proton and net proton multiplicities 

Using the pr integrated density distributions dn/dxp extracted above for the inclusive 
p-i-C interaction as well as for the different sub-components, the corresponding baryonic multi- 
plicities may now be determined by integration over Xi?. This results in the following numbers: 



inclusive proton multiplicity (rip) = 1.881 (29) 

inclusive anti-proton multiplicity (rip) = 0.052 (30) 

multiplicity of pair produced protons (rip) = 0.0653 (31) 

inclusive net proton multiplicity (n^^^) = 1.816 (32) 

There are two predictions respectively for the net projectile component and for the ex- 
pected net target component established from the elementary p-i-p collisions [i3J: 

predicted net projectile component (72^''°^'^'''^'') = 0.550 (33) 

This assumes baryon number conservation and a negligible increase in neutron and hy- 
peron production. 

predicted net target component ^^targ,pred^^ ^ g ggg ^^4) 

This uses isospin invariance and the measured number of 1.6 projectile collisions. 
Subtracting the predicted target component from the total net proton density in the pro- 
jectile fragmentation region, the multiplicity of net projectile protons is measured as 

measured net projectile component (raP''"'''") = 0.520 (35) 

which reproduces the expected multiplicity to within 5.5% and leaves room for some increase 
of the hyperon yields in p-i-C collisions. 
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Subtracting the predicted target component from the total net proton density in the back- 
ward region, the nuclear component is measured as 

measured nuclear component (n™^'''"'^^-^) = 0.655 (36) 

A Gauss fit to the nuclear component in the symmetric regions of the target-subtracted 
proton density yields the multiplicity 

Gaussian nuclear component = 0.631 (37) 

which complies to within 3.6% with the straight-forward subtraction. 

Finally there is the modified target component obtained by subtraction of the Gaussian 
nuclear component from the overall proton density which smooths the diffractive peak in accor- 
dance with the Fermi motion of the hit nucleons. This modified target component results in the 
multiplicity 

modified target component (^'^g'™'^) = 0.670 (38) 

which reproduces the predicted multiplicity (l34l) to within 2.9%. 

In conclusion it may be stated that the total measured net proton multiplicity of 1.82 per 
inelastic event is made up by the superposition of three components of comparable magnitude, 
namely a projectile multiplicity of 0.52, a target multiplicity of 0.67 and a nuclear component 
of 0.63 units. The multiplicities of anti-protons and pair produced protons are 0.052 and 0.0653, 
respectively, corresponding to 2.9% and 3.6% of the total net proton yield. 

14 Proton production in p+C collisions: pr dependence 

Following the study of px integrated yields in the preceding section, the double- 
differential proton cross sections will now be studied for the three components of projectile 
fragmentation, target fragmentation, and nuclear cascading, thus adding the transverse dimen- 
sion to the experimental scrutiny. In order to comply with the discussion of the integrated proton 
yields dn/dxp. Sect. \T3\ the double differential cross sections will be used in their non-invariant 
form 

J__(Pn_ ^ 2n ^f{xF,PT) ^29) 

PrdXpdpT (Tinel 2 E 

see also (flOl) . In a first step, the net proton cross sections will be defined using the results on 
anti-protons. Sect. [121 A detailed comparison of the overall backward cross sections with the 
predicted contribution from target fragmentation will reveal a necessary px dependent modi- 
fication of this component which will then be employed to extract in turn the projectile and 
nuclear components in their px dependence. 

14.1 Pair produced protons and net proton density 

Using the results of the discussion of the double differential anti-proton yields and their 
separation into target and projectile components in Sect. 112. 21 the yield of pair produced protons 
may be estimated by multiplying the projectile component by the isospin factor 1 .6, maintaining 
the Px enhancement as parametrized in (l25l) . The target component stays equal to the anti-proton 
yield due to isospin symmetry. The total and the resulting net proton densities are shown in 
Fig. ll02l as a function of for different values of px- 
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Figure 102: Total and net proton density — dxldpT ^ function of xp for different values 
indicated in the panel. The distribution at pt = 0.2 GeV/c is multiplied by factor of 1 .2 for better 
separation 



14.2 Target component 

As discussed in Sect. 113. 2^ the target component of net proton production in p+C inter- 
actions may be predicted from p+p collisions using the number of projectile collisions (u) = 1 .6, 
the target overlap function R'^'^^^ (l22l) and Fig. |7T1 and the measured proton and neutron cross 
sections [3|. The only additional assumption to be made for the double differential yield is the 
equality of the pr distributions for neutrons and protons, see [|3]| for experimental evidence. The 
result for the target component and its subtraction from the overall proton yield is shown in 
Fig. [Ml 

Compared to the same prediction and subtraction for the px integrated densities, Figs.|94l 
and|95l a pt dependence becomes visible. It is the region of xp between -0.2 and -0.3 which 
is extremely sensitive to the predicted target yield. At transverse momenta below the mean pt 
(Fig.[55l) the prediction overshoots the total density by about 20%, whereas for higher it falls 
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Figure 103: Double differential net proton yield — , as a function of xp for five different 
values of pt- Full line: total yield, broken line: predicted target component, and dotted line: 
after subtraction of the target component 



low by up to 40% at = 1 GeV/c. In this xp region neither the possible projectile feed-over 
nor the nuclear component may explain, by their limited x p range, the observed pp dependence. 
Instead a yield suppression at low pp followed by an enhancement at high pp similar to the one 
found for anti-protons. Sect. 112.21 has to be invoked also for the target contribution. In order to 
preserve a physical dependence of the projectile and nuclear components, see Sects. [143] and 
1 14.41 below, an xp dependent pattern of the target enhancement as shown in Fig. 11041 imposes 
itself. 




Figure 104: Enhancement factor of the target component a) as a function of for pp values 
between 0.2 and 1 GeV/c and b) as a function of pp for xp = -0.25. The shaded region represents 
the error margins 
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This phenomenon is also borne out by a comparison of the mean pt values between p+C 
and p+p interactions shown in Fig. 11051 




Figure 105: Comparison of the mean pt values of protons in p+C (open circles) and p+p inter- 
actions (closed circles) as a function ofxp from xp = -'^■Q xp = +0.95 



If the target component would show a px dependence equivalent to p+p interactions, 
the mean px values should approach the ones for p+p dX xp around -0.2 to -0.3 where target 
fragmentation prevails. 

Applying this enhancement, the target component and its subtraction from the total mea- 
sured yield take the pattern shown in Fig. 11061 
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Figure 106: Double differential net proton yield — , as a function of xp for five values of 
Pt- Full line: total yield, broken line: predicted enhanced target component, dotted line: after 
subtraction of the enhanced target component 
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This subtraction leaves clearly defined and separated projectile and nuclear components 
which will be discussed in the subsequent sections 



14.3 The projectile component 

The hadronization of the projectile is already clearly visible in the dotted line of Fig. [T06l 
It saturates the total yield at xp > 0.2 and comes down to zero at x^? ~ -0.2 due to the limited 
range of the baryonic overlap function. In order to put this behaviour in perspective as far as 
the X p distributions for different transverse momenta are concerned, the projectile components 
normalized to their maximum densities are presented in Fig. 11071 




Figure 107: Net proton density — ^^'^^^^ from projectile fragmentation as a function of xp 
normalized to its maximum value for five px values between 0.2 and 1 GeV/c 



Compared to the px integrated distribution, Fig.|95l there is a clear dependence of baryon 
number transfer on px- The maximum density shifts from xp = 0.3 at low p^- to s^r = 0.1 at 
1 GeV/c. The density atxp> 0.6 corresponding to low mass excitation is strongly reduced with 
increasing px whereas the transfer function is successively extending further into the backward 
hemisphere. 

Further information concerning the px dependence comes from a direct comparison to 
p+p interactions. Extracting the projectile component from the symmetric p+p collision using 
the baryonic overlap function determined in Sect. II 1.1] the ratio 

/I d^n / f 1 <fn y 

\pxdxpdpx) I yppdxpdpx) ^ ^ 

may be obtained as shown in Fig. llOSI as a function of px for different values of xp. 

A clear increase with px of the net proton yields is visible for all values of xp. The 
overall increase of the cross sections with decreasing a;i? is due to the general baryon number 
transfer in p+C interactions visible in the px integrated distributions of Fig. [57] ("minimum bias" 
compared to p+p). This can be normalized out by multiplying by the integrated density 
ratio between p+p and p+C collisions: 
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Figure 108: Density ratio RF^°^ of the projectile components of p+C and p+p interactions as a 
function of px for different values of xp 



The resulting ratio -Rn^m shown in Fig. 11091 again as a function of px for different 
values of xp- 




^■^0 0.5 1 1.5 

[GeV/c] 

Figure 109: Normalized density ratio R^"j^ as a function of px for different values of Xi? 
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The wide spread observed in Fig. 11081 reduces to a narrow band of pr enhancements 
which are about independent on xp except for the xp range approaching diffraction. Again, as 
for p in Sect. 1 12.2] and for the target component in the preceding section, the normalized densi- 
ties are reduced by 10-20% at low px and enhanced to factors of about 1.5 at = 1.5 GeV/c. 



14.4 The nuclear component 



The onset of the contribution from nuclear cascading is already visible in Fig. 1 1061 by 
the dotted lines at xp < -0.25. This range is extended to the far backward hemisphere down to 
Xp = -1.6 in Fig. II 101 In this Figure the total proton density, the predicted target contribution 
and the target subtracted yield are shown as a function of for five values of pp between 0.2 
and 1 GeV/c. 
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Figure 1 10: Double differential net proton yield as a function ofxp for five values of pp. Broken 
line: total proton yield, dotted line: predicted target contribution, full line: target subtracted 
proton density 



As already discussed for the pp integrated densities in Sect. 113.41 (Figs. |97] and l98l). the 
presence of the diffractive peak from proton fragmentation in the predicted target contribution 
creates a sharp spike between xp ~ -0.9 and xp = -\.0 which should be visible in the total proton 
density distribution in this area. The absence of such structure in the measured yield indicates, 
as argued in Sect. 113.41 a smoothing of the predicted structure via quasi-elastic re-scattering of 
the diffrative protons inside the Carbon nucleus. This re-scattering extends typically up to ±0.2 
units ofxp around xp = -l [|45l . Beyond this range the target subtracted yield shows a Gaussian 
behaviour as presented in Fig. II 1 II 

The fit parameters show, within tight errors, a stable mean value in xp between -0.92 
and -0.93 in agreement with the pp integrated fit (l28l) . The rms deviation increases with pp from 
0.22 to 0.3 units of xp. The pp dependences of the maximum density and of the yield integrated 
over X p are shown in Fig. 11121 
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Figure 111: Double differential, target subtracted proton yields as a function of xp for -1.6 < 
xp < -0.3 with the exclusion of the range -1.1 < xp < -0.8 for five values of pr between 
0.2 and 1.0 GeV/c. Superimposed as full lines are Gaussian fits with mean values (xp) and a 
indicated in each panel. The last panel shows plots of {xp) (left scale) and a (right scale) as a 
function of pt 




[GeV/c] [GeV/c] 

Figure 112: a) Maximum proton density and b) yield integrated over of the nuclear 

component as a function of pr- The full lines represent data interpolations, the broken lines 
correspond to the pt dependence of the proton density in p-i-p interactions at \xp\ = 0.95 [|3||, 
normalized at = 0.4 GeV/c 



It is interesting to observe that the pt dependence of the nuclear component is not Gaus- 
sian and corresponds to the one measured in the diffractive region of p-i-p interactions jSl. 

Using the Gaussian fits as an estimator of the nuclear component a modified target com- 
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ponent may be obtained by subtracting these fitted densities from the total proton yield as shown 
in Fig. [ml 




L5 



Figure 113: Double differential proton yields as a function of xp for five values of pt- Broken 
lines: total proton yield, dotted lines: Gaussian fit to the nuclear component, full lines: modified 
target component obtained by the subtraction of the Gaussian fits from the total proton density 

The resulting modified target component extends now toxp values below -1.0 and shows 
a smooth behaviour through the region of proton diffraction in accordance with the px integrated 
yield, Sect.[TT4l 

15 Conclusion 

As part of a comprehensive study of hadronic production in elementary and nuclear col- 
lisions at the SPS, new data from the NA49 experiment on proton, anti-proton, neutron and light 
ion production in minimum bias p-i-C interactions are presented. Making full use of the accep- 
tance coverage and the particle identification capabilities of the NA49 detector, a wide phase 
space area from the far forward direction at Feynman of 0.95 to the far backward direction 
down to Xi;' = -0.8 has been exploited. Using available data from a Fermilab experiment the data 
coverage could be further extended into the nuclear fragmentation region down to Xi? ~ -2.0. 
In addition, deuteron and triton production have been studied in the lab momentum range from 
0.25 to 3 GeV/c making available there for the first time cross sections in the low to medium px 
region. In addition a limited amount of data with grey proton detection allows for a first look at 
the centrality dependence of baryon number transfer. 

Given the complete phase space coverage of the combined data set, the main aim of this 
publication is the separation and isolation of the three components of hadronization in p+C col- 
lisions, namely projectile fragmentation, target fragmentation and nuclear cascading. This study 
has been conducted both for integrated quantities and using double differential cross sections 
to obtain a complete view of the px dependence. For this aim, the baryonic overlap functions 
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from the projectile to the target regions and vice versa have been determined experimentally 
using both the elementary p+p and the asymmetric p+C interactions. This has been achieved in 
a completely model-independent way relying essentially on baryon number conservation and 
isospin symmetry. For anti-proton production the absence of a nuclear contribution has been 
shown and the superposition of the target and the projectile fragmentation has been established 
using the known number of projectile collisions inside the Carbon nucleus. This allows for the 
definition of the yield of pair produced protons and thereby of net proton densities. 

In contrast to the a priory unknown projectile and nuclear components, the target frag- 
mentation occupies a special place as it should at least to first order be describable by a super- 
position of single nucleon fragmentations taking of course into account the number of projectile 
collisions and isospin effects. This approach works out well for pt integrated proton densities 
with the exception of the diffractive contribution contained in the elementary interactions. This 
contribution is evidently smeared out by secondary, quasi-elastic interactions of the correspond- 
ing low momentum protons in the nuclear rest system. The subtraction of the thus predicted 
target component yields both the projectile and the nuclear components. The former is proven 
to preserve baryon number combined with a sizeable amount of baryon transfer of order 0.15 
units of xf- The latter turns out to have a Gaussian shape in centred at = -0.92 with a 
substantial FWHM of 0.6 units far in excess of the narrow xp distribution centred close to -1.0 
which might be expected from Fermi motion alone. In the case of the p-i-C collisions studied 
here, the three components have comparable total yields of 0.52, 0.67 and 0.65 net protons re- 
spectively for the projectile, target and nuclear contributions. This sums up to 1.84 net protons 
in total which complies with the direct yield integration to within 1%. 

The study of double differential cross sections gives access to the additional dimension 
of transverse momentum. All three components show a distinct transverse activity which goes 
beyond the naive expectation from elementary collisions or nuclear binding. The projectile com- 
ponent features, both for protons and anti-protons, a suppression by 10-20% at low pt followed 
by an enhancement of about 50% at pt= 1-5 GeV/c. This pattern is rather a; i? independent in the 
forward hemisphere. The target component shows a similar behaviour, however with a distinct 
xp dependence centred dXxp ^ -0.3. The nuclear component finally has a wide, non-Gaussian 
Pt dependence which goes far beyond the one expected from Fermi motion and which is shown 
to be comparable to the one measured in the diffractive region of p+p interactions. 

As for the previous publications the measured cross sections and two-dimensional 
interpolation values are available on the Web Site ll46l . 
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